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Foreword 

This Ethiopian Standard has been prepared under the direction of the Technical Committee for Building 
structures and elements of building (TC 41) and published by the Ethiopian Standards Agency (ESA). 

The standard is identical with ISO 15099:2003 Thermal performance of windows, doors and shading 
devices — Detailed calculations published by International organization for Standards (ISO). 

For the purpose of this Ethiopian Standard , the adopted text shall be modified as follows. 

• The phrase "International Standard" shall be read as "Ethiopian Standard"; and 

• A full stop (.) shall substitute comma (,) as decimal marker. 
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Thermal performance of windows, doors and shading 
devices — Detailed calculations 



1 Scope 

This International Standard specifies detailed calculation procedures for determining the thermal and optical 
transmission properties (e.g., thermal transmittance, total solar energy transmittance) of window and door 
systems based on the most up-to-date algorithms and methods, and the relevant solar and thermal properties 
of all components. 

Products covered by this International Standard include windows and doors incorporating: 

a) single and multiple glazed fenestration products with or without solar reflective, low-emissivity coatings 
and suspended plastic films; 

b) glazing systems with pane spacing of any width containing gases or mixtures of gases; 

c) metallic or non-metallic spacers; 

d) frames of any material and design; 

e) fenestration products tilted at any angle; 

f) shading devices; 

g) projecting products. 

2 Normative references 

The following referenced documents are indispensable for the application of this document. For dated 
references, only the edition cited applies. For undated references, the latest edition of the referenced 
document (including any amendments) applies. 

ISO 7345, Thermal insulation — Physical quantities and definitions 

ISO 8301, Thermal insulation — Determination of steady-state thermal resistance and related properties — 
Heat flow meter apparatus 

ISO 8302, Thermal insulation — Determination of steady-state thermal resistance and related properties — 
Guarded hot plate apparatus 

ISO 9050, Glass in building — Determination of light transmittance, solar direct transmittance, total solar 
energy transmittance, ultraviolet transmittance and related glazing factors 

ISO 9288, Thermal insulation — Heat transfer by radiation — Physical quantities and definitions 

ISO 9845-1, Solar energy — Reference solar spectral irradiance at the ground at different receiving 
conditions — Part 1: Direct normal and hemispherical solar irradiance for air mass 1,5 
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ISO 10077-2:2003, Thermal performance of windows, doors and shutters — Calculation of thermal 
transmittance — Part 2: Numerical method for frames 

ISO 10211-1, Thermal bridges in building construction — Heat flows and surface temperatures, Parti: 
General calculation methods 

ISO/CIE 10526:1999, CIE standard llluminants for colorimetry 

ISO/CIE 10527, CIE standard colorimetric observers 

ISO 12567-1 , Thermal performance of windows and doors — Determination of thermal transmittance by hot 
box method — Part 1: Complete windows and doors 

EN 12898, Glass in building — Determination of the emissivity 



3 Symbols 



3.1 General 



Symbols and units used are in accordance with ISO 7345 and ISO 9288. The terms, which are specific to this 
International Standard, are listed in Table 1. 



3.2 Symbols and units 



Table 1 — Terms with their symbols and units 



Symbol 


Term 


Unit 


A 


area 


m 2 


A i 


portion of absorbed solar energy by the rth glazing layer 


1 


A R 


aspect ratio 


1 


b 


width (breadth) of a groove or slit 


mm 


c p 


specific heat capacity at constant pressure 


J/(kgK) 


d 


thickness 


m 


d , 


thickness of glazing cavity 


m 


E 


irradiance 


W/m 2 


E S (A) 


solar spectral irradiance function (see ISO 9845-1) 


1 


E y {X) 


colorimetric illuminance (CIE D65 function in ISO/CIE 10526:1999) 


Ix 


g 


acceleration due to gravity 


m/s 2 


G 


parameter used in the calculation of convective heat transfer coefficients; see Equation (48) 


1 


h 


surface coefficient of heat transfer 


W/(m 2 K) 


H 


height of glazing cavity 


m 


I 


total density of heat flow rate of incident solar radiation 


W/m 2 


/"(A) 


spectral heat flow rate of radiant solar energy between rth and z'+1th glazing layers 
travelling in the external ( + ) or internal (~) direction 


W 
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Table 1 (continued) 



Symbol 


Term 


Unit 


J 


radiosity 


W/m 2 


I 


length 


m 


M 


molecular mass 


mole 


N 


number of glazings + 2 


1 


Nu 


Nusselt number 


1 


P 


pressure 


Pa 


1 


density of heat flow rate 


W/m 2 


r 


reflectance: portion of incident radiation reflected such that the angle of reflection is equal 
to the angle of incidence 


1 


R 


thermal resistance 


m 2 K/W 


R(X) 


photopic response of the eye (see ISO/CIE 10527) 




iff 


universal gas constant 


J/(kmolK) 


Ra 


Rayleigh number 


1 


Ra x 


Rayleigh number based on length dimension x 


1 


Si 


density of heat flow rate of absorbed solar radiation at z'th glazing laver 


W/m 2 


f perp 


largest dimension of frame cavity perpendicular to heat flow 


m 


T 


thermodynamic temperature 


K 


AT i 


temperature drop across z'th glazing cavity, AT; = \T U - r b/+1 | 


K 


u 


air velocity near a surface 


m/s 


U 


thermal transmittance 


W/(m 2 K) 


V 


free-stream air speed near window, mean air velocity in a gap 


m/s 


x,y 


dimensions in a Cartesian co-ordinate system 


1 


z 


pressure loss factor 


1 


a 


absorption 


1 


P 


thermal expansion coefficient of fill gas 


K- 1 


£ 


total hemispherical emissivity 


1 


r 


angle 


o 


9 


temperature 


°C 


a 


Stefan-Boltzmann constant, 5,669 3 x 1CT 8 


W/(m 2 -K 4 ) 


X 


thermal conductivity 


W/(mK) 


K 


wavelength 


m 


v 


dynamic viscosity 


Pas 


p 


density 


kg/m 3 


T 


transmittance 


1 


T s 


total solar energy transmittance: the portion of radiant solar energy incident on the 
projected area of a fenestration product or component that becomes heat gain in the 
internal conditioned space 


1 


4> 


parameter used in the calculation of viscosity and of thermal conductivity; 
see Equations (62) and (67) 


1 


<P 


function used in the calculation of heat transfer; see Equation (112) 


1 





heat flow rate 


W 


w 


linear thermal transmittance 


W/(mK) 
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3.3 Subscripts 

The subscripts given in Table 2 shall be applied. 



Table 2 — 


Subscripts and meanings 


Subscript 


Meaning 


ai 


air 


av 


average 


b 


backward 


bo 


bottom of a gap 


cc 


condition on the cold side 


cdv 


conduction/convection (unvented) 


eg 


centre of glass 


ch 


condition on the hot (warm) side 


cr 


critical 


cv 


convection 


de 


divider edge glass 


dif 


diffuse 


dir 


direct 


div 


divider 


eff 


effective 


eg 


edge of glass 


eq 


equivalent 


ex 


external 


f 


frame 


fr 


frame (using the alternative approach) 


ft 


front 


gv 


glass or vision portion 


ht 


hot 


hz 


horizontal 


i 


counter 


int 


internal 


inl 


inlet of a gap 


J 


counter 


m 


mean 


mix 


mixture 


n 


counter 


ne 


environmental (external) 


ni 


environmental (internal) 


out 


outlet of a gap 


P 


panel 


r 


radiation or radiant 


red 


reduced radiation 


s 


surface 


sc 


source 


sk 


sink 


si 


solar 


t 


total 


tp 


top of a gap 


V 


number of gases in a gas mixture 


V 


vertical 


z 


at distance z 


¥ 


perimeter 


2D 


coupling 
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4 Determination of total window and door system properties 
4.1 Thermal transmittance 

4.1.1 General 

This International Standard presents procedures by which detailed computations can be used to determine 
the thermal transmission properties of various product components, which are then used to determine the 
thermal transmission properties of the total product. Where national standards allow, test procedures may be 
used to determine component and total product properties. 

The total properties for window and door products are calculated by combining the various component 
properties weighted by either their respective projected areas or visible perimeter. The total properties are 
each based on total projected area occupied by the product, A t . The projected component areas and the 
visible perimeter are shown in Figure 1 . 




r 

> 

en 



Key 

1 perimeter length at sight line 



Figure 1 — Schematic diagram showing the window projected areas and vision perimeter 

Clause 4 describes the procedure for calculating thermal transmittance, total solar transmittance and visible 
transmittance for the complete product. 4.1 describes the procedure for calculating thermal transmittance. The 
effect of three-dimensional heat transfer in frames and glazing units is not considered. 4.1.4 describes an 
alternative procedure for calculating edge of glass and frame thermal indices U de , U e „, U t and E/ fr , which are 
used in area-based calculations. Clause 5 describes the procedure for calculating the required centre-glass 
properties r sgv and r gv . Clause 6 describes the procedure for calculating the linear thermal transmittance, f, 
which accounts for the interaction between frame and glazing or opaque panel. Clause 7 contains the 
procedure for dealing with shading devices and ventilated windows. Clause 8 describes the procedure for 
determining and applying boundary conditions. The thermal transmittance of the fenestration product is given 
by: 



U, 



I.Ag V U gy +I.A 1 U 1 +Z/yiP 



(1) 



where A gv and A f are the projected vision area and frame area, respectively. The length of the vision area 
perimeter is Aj,, and ^is a linear thermal transmittance that accounts for the interaction between frame and 
glazing or the interaction between frame and opaque panel (e.g., a spandrel panel). 
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The summations included in Equation (1) are used to account for the various sections of one particular 
component type; e.g. several values of A^ are needed to sum the contributions of different values of U f 
corresponding to sill, head, dividers and side jambs. 

Figure 2 illustrates the division into components for the alternative approach described in 4.1.4, in which the 
edge-of-glass and divider-edge areas are 63,5 mm (2,5 in) wide. The sum of all component areas equals the 
total projected fenestration product area. 
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Figure 2 — Centre-of-glass, edge-of-glass, divider, divider-edge, 
and frame areas for a typical fenestration product 

4.1.2 Glazed area thermal transmittance 

The thermal transmittance can be found by simulating a single environmental condition involving 
internal/external temperature difference, with or without incident solar radiation. Without solar radiation, the 
thermal transmittance is the reciprocal of the total thermal resistance. 



U 



gv 






(2) 



and when solar radiation is considered, then: 



U 



9int(/ s =0) 



gv 



(3) 
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where g int (/ s = 0) is the net density of heat flow rate through the window or door system to the internal 
environment for the specified conditions, but without incident solar radiation, in W/m 2 . The condition "without 
solar radiation" is used because all effects on the thermal resistances due to incident solar radiation are 
incorporated in the total solar energy transmittance or r s -value [see Equation (14)], and f ni and T ne are the 
environmental temperatures, as defined in Equation (7). 

R t is found by summing the thermal resistances at the external and internal boundaries, and thermal the 
resistances of glazing cavities and glazing layers. See Figure 3. 



* n n 

*t= — + 2>+2>gv,,- 



1 



'int 



(4) 



"ex i=2 i=\ 

where the thermal resistance of the /th glazing is: 



R 



gv,r 



gv,r 



A 



(5) 



gv,*' 



and the thermal resistance of the /th space, where the first space is external environment, the last space is 
internal environment and the spaces in between are glazing cavities, (see Figure 3): 



R< 



1 f,i - 1 W-1 



(6) 



where r f ,-, and T b M are the external and internal facing surface temperature of the /th glazing layer. 

The environmental temperature [as defined in Equation (7)] is a weighted average of the ambient air 
temperature and the mean radiant temperature, T m , which is determined for external and internal 
environment boundary conditions (see boundary conditions in 8.4.1). 




1 
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/ + 1 



c, int 



n + 1 
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int 



rm, int 



Key 

1 gap 

2 glazing 

Figure 3 — Numbering system for glazing system layers 

The environmental temperature, T n , is: 

T _ "cv 1 ai T "r 1 rm 
n ~ , , 

«cv +«r 

where h cv and h r are determined according to the procedure given in Clause 8. 



(7) 
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4.1.3 Frame area/edge-glass thermal indices 

In order to convert the results of a two-dimensional numerical analysis to thermal transmittances, it is 
necessary to record the rate of heat transfer from the internal environment to the frame and edge-glass 
surfaces (in the absence of solar radiation). The linear thermal transmittance, W, values and frame thermal 
transmittances shall be calculated according to the following equations. 

<? = L 2D -U f l f -U gv l gv (8) 

where L 2D is thermal coupling coefficient determined from the actual fenestration system. 

Ui = L p 2D -U p l p (g) 

h 

where 

L p 2D is thermal coupling coefficient determined from the frame/panel insert system; 

Up is the thermal transmittance of foam insert; 

L is the internal side exposed length of foam insert (minimum 100 mm); 

l f is the internal side projected length of the frame section; 

/ gv is the internal side projected length of the glass section (see Figures C.1 and C.2 of 
ISO 1 0077-2:2003, for further details on the definition of / fr and / p ). 

The detailed procedure for determining L 2D is also given in ISO 10211-1. 

4.1.4 Alternative approach (see Figure 2) 

An alternative method is available for calculating frame thermal transmittance, U fr Using this method it is 
unnecessary to determine the linear thermal transmittance, X F. Instead, the glass area, A gy , is divided into 
centre-glass area, A c , plus edge-glass area, A e , and one additional thermal transmittance, C/ eg , is used to 
characterize the edge-glass area. If dividers are present then divider area, ^ div and divider thermal 
transmittance, U dw are calculated, as well as corresponding divider edge area, A de and thermal transmittance, 
C/ de . The following equation shall be used to calculate the total thermal transmittance: 

TT X^cg^c+X^fr^f + X (7 eg^e+X ;7 div^div+X ;7 de^de .... 

U t = (1U) 

A \ 
where t/ fr , and U eg can be determined from the following equations: 

t/ fr = 0fr r (11) 

h (Tni - ^ne ) 

U eQ = 7—^ v ( 12 ) 

e " 7 It . -T \ v ' 

'egl/ni ^nej 

and where / f is projected length of frame area and / eg is the length of edge of glass area and is equal to 
63,5 mm. These lengths are measured on the internaf side. The quantities c/> fr and cP eg are heat flow rates 
through frame and edge-glass areas (internal surfaces), respectively, including the effect of glass and spacer, 
and both are expressed per length of frame or edge-glass. The calculations shall be performed for each 
combination of frame and glazing with different spacer bars. 
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The summations included in Equation (10) are used to account for the various sections of one particular 
component type; e.g., several values of A f must be used to sum the contributions of different values of U fr 
corresponding to sill, head and side jambs. 

It should be noted that the two different approaches entail different definitions of frame thermal transmittance, 
denoted £/ f and U fr The primary difference is that the (7 fr includes the some of the heat transfer caused by the 
edge seal, whereas (7 f does not. The comparison of frame properties for two different products is only 
meaningful if the same calculation procedure has been used in both cases. 

The U { values for windows calculated by the two methods may differ because of differences in the way frame 
and edge heat transfer is treated at the corners, particularly because the three dimensional effects are 
neglected. This difference is more pronounced for smaller windows. The choice of / eg = 63,5 mm is made to 
reduce the discrepancy between the two alternative approaches. 

4.2 Total solar energy transmittance 

4.2.1 General 

The total solar energy transmittance of the total fenestration product is: 



A t 



(13) 



where t q and r f are the individual total solar energy transmittance values of the vision area and frame area, 
respectively. The summations are included for the same reason that they appear in Equation (1) and shall be 
applied in the same manner to account for differing sections of one particular component type. 

NOTE Equation (13) includes an assumption that the solar transmittance of the edge of glass is the same as that of 

the centre of glass area. 

4.2.2 Vision area total solar energy transmittance 

The total solar energy transmittance can be determined for conditions involving internal/external temperature 
difference and any level of incident solar radiation. It is found by calculating the difference between the net 
heat flow rate into the internal environment with and without incident solar radiation. 

(14) 



2"S = 


q 


int ' 


-9int( 7 s =0) 






's 


where 








?int 






is the net 



is the net density of heat flow rate through the window or door system to the internal 
environment for the specified conditions, in W/m 2 ; 

g int (/ s = 0) is the net density of heat flow rate through the window or door system to the internal 
environment for the specified conditions, but without incident solar radiation, in W/m 2 . 

For the equivalent expression for U, see Equation (3). 

The net density of heat flow rates, q mt and q mi (I s = 0) are calculated in 5.3.1 [Equation (27), for index i = int]. 

For a glazing assembly in which a shading device is involved, the amendments to the equations of 5.2 as 
given in 7.2 shall be applied. 
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4.2.3 Frame total solar energy transmittance 

The frame total solar energy transmittance shall be calculated using the approximate equation: 



r f =« f 



Ui 



A s 



(15) 



where A s is the developed surface area. 

The external surface heat transfer coefficient (combined convective/radiative) at the frame, h ex , is 
"ex = "cv.ex + "r,ex- 

If the alternative method of calculating U t is being used, U fr should be used instead of (7 f in Equation (15). 

More detailed two-dimensional or three-dimensional calculations, including the effects of off-normal solar 
radiation, shading, reflected solar radiation and solar radiation transmitted to the internal frame surfaces, can 
be performed in a manner analogous to Equation (14), and subject to boundary conditions given in 8.6. 

4.3 Visible transmittance 

The visible transmittance of the total fenestration product is: 



A, 



(16) 



5 Vision area properties 
5.1 Glazing layer optics 

5.1.1 General 

For glazing units only, the optical properties can be determined using ISO 9050. Clause 7 contains the 
extensions needed to model vented windows. 

5.1.2 Solar 

The solar optical properties needed to describe the z'th glazing are: a) the front (external side) spectral 
reflectance, r ftl (/l w ); b) the back (internal side) spectral reflectance, r b ,j(/l w ); and c) the spectral 
transmittance, z;'(/l w ). See Figure 4. 

NOTE More information about r ft . (A w ), r b .(/L w ) and tj(A^) can be found in [2]. 




















t,(X) 
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1 outdoor side 

2 z'th glazing layer 

3 indoor side 



Figure 4 — Outdoor and internal spectral transmittance of a glazing layer 
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The solar optical data shall be measured in accordance with ISO 9050. Intermediate values of r u (A^), 
or t^X^) are found by linear interpolation. 



,(4J 



5.1.3 Long-wave 

The long-wave optical properties needed to describe the /th glazing are: a) the front (external side) 
hemispheric emissivity, s^{, b) the back (internal side) hemispheric emissivity, s bi ; and c) the hemispheric- 
hemispheric transmittance,' r t . These total optical properties apply to wavelengths from 5 urn to 50 urn. 

The long-wave reflectance data shall be measured in accordance with EN 12898. Values of the normal 
emissivity resulting from this procedure shall be converted to hemispherical emissivity using the procedure 
described in [3] or in EN 12898. The integration needed to convert measured spectral data to the required 
total longwave optical properties, £- ft t , s b t and r t shall be carried out in accordance with [3] or EN 12898. 

Some windows are constructed with suspended or stretched layers of thin plastic film between glass panes to 
make triple or quadruple glazing. When these layers are covered with a low-em issivity coating they are 
generally opaque in the infrared, so that r t = and hemispherical emissivity can be calculated as in [3]. For 
partially transparent films such as polyethylene terephthalate (PET), both specular transmittance and 
reflectance should be measured. Using the bulk model, the optical indices of the material can then be 
calculated and used to derive the hemispherical properties. 



NOTE 



See [4] for more information. 



5.2 Glazing system optics 

5.2.1 Spectral quantities 

The path of incident solar radiation within the various layers of the glazing system shall be modelled by the 
methods described in ISO 9050 or by any other exact method. 

NOTE Depending on future modifications of ISO 9050, specific additions may be added to this International Standard 

covering the effects of optical properties of products (shading devices, diffusing panes, etc.) not adequately covered by 
ISO 9050. 




a, {X) 



M*) 
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Figure 5 — Absorption of the /th glazing layer and solar spectral transmittance 
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Figure 5 shows how a window with n glazing layers together with the external = 0) and internal (i = n+ 1) 
spaces can be treated as an n + 2 element array. It is necessary to determine the portion of incident solar 
radiation, at a given wavelength, that is absorbed at each of the glazing layers. This quantity is denoted a t {X w ) 
at the z'th glazing layer. Similarly, it is necessary to determine the solar spectral transmittance of the glazing 
system, t s {XJ. 

These quantities, a,(/l w ) and r s {X w ), shall be calculated in accordance with ISO 9050 while setting the 
reflectance of the conditioned space to zero. Any other method that can be shown to provide an exact solution 
is acceptable. 

NOTE The solution technique described in [5] is summarized in Annex A. 

5.2.2 The solar spectrum 

The spectral distribution of the incident solar radiation, E(X W ), is needed to calculate total optical properties 
and various total energy flow. Values of E(X VJ ) are reported at N s \ values of X^ (denoted here as E(X W ) and 
Xfij, respectively). Intermediate values of E(X W ) shall be found by linear interpolation of the tabulated values. 

5.2.3 Absorbed amounts of solar radiation 

The total flow rate of solar radiation absorbed at the /th glazing layer, S t , is determined by numerical 
integration over the solar spectrum according to Equations (16), (17) and (18). 



JVsl-1 
X a i[ X ™ j I y+ij^sl^w j/ y+l) A/l w/ 

Ai=J ^^ 07) 

X ^Sl^Wj/j+lj^Wj 
7=1 

AX W j = X W j + -\ - X W j (18) 

S i =A t xI a (19) 

where A i is the portion of the total incident solar radiation (on the glazing system) that is absorbed by the z'th 
glazing layer, and a t (^,/, + i) is the value of a { that is representative of the wavelength band from X^j to /l w+1 
and is given by 

a «'( lw iO'+i) = I a, '(' lw i) + /''( /lw i+i) ( 20 ) 

and 

E *\*vjlj+V- 2 { ' 

Values of E^X^) are given in ISO 9845-1. 
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5.2.4 Solar transmittance 

The solar transmittance of the glazing system is: 



^Sl = Jf— A/l W j = ^W y+1 ~ ^W y (22) 

X ^(^Wy/y+lj^Wy 

where ^{'^yy+i) is given by Equation (21) and 

T S\[ ;l Wjlj + -\)= 2 r Sl( /l Wy)+2 r Sl(' ;l Wy + l) (23) 

5.2.5 Visible transmittance 

Visible transmittance, r vs , is calculated using a weighting function that represents the photopic response of the 
eye, ^(/LJ. R{A w j} is tabulated for N vs values of A^. r vs is given by: 



w vs-1 
7=1 



AA W , - A w , +1 - A w , (24) 



X ^VS^Wy/y+lj^^Wy/y+lJ^Wy 

7=1 



where 



^^wy/y + ij- 2 (25) 

^VsKy/y+^^vs^Wyj + ^vs^Wy+l) (26) 

Values of £ v (/l w ) are given in ISO/CIE 10526. 
and ^ S |(^ w/ 7 /+ i) is given by Equation (23). 

5.3 Vision area heat transfer 

5.3.1 Glazing layer energy balance 

Longwave radiative exchange between glazing layers and conductive heat transfer within each glazing layer 
can be described using fundamental relations. Calculations dealing with convective heat transfer depend upon 
correlations based on experimental data. 
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Figure 6 — Energy balance on glazing layer i 



Figure 6 shows the /th glazing in a sloped multilayer array. The values of four variables are sought at each 
glazing. These are the temperatures of the external and internal facing surfaces, T fti and T bj plus the radiant 
heat leaving the front and back facing surfaces (i.e., the radiosities), J ni and J bj . In' terms of these variables, 

9i is: 



ii = h c\i,i [ T n,i ~ Tb.j-i ) + J n,i ~ J b,i-\ 
The solution is generated by applying the following four equations at each glazing: 

J n,i = £ n,i (7T n,i + T i J fu+-\ + r ft,i J b,i--\ 

J b,i = s b,i aT b,i + T i J b,i^ + r b,i J ft,i+-\ 



'gv,i 



2/1 



(2q M +Si) 



(27) 

(28) 
(29) 
(30) 

(31) 



gv,i 



A by-product of the analysis is the temperature profile through each glazing. 



*/(*)< 



Si 



^-^gv.i'^gv.r 



'fU 



\ 



tc 



2l c 



z + T h 



(32) 



'gv,; ^- /L gv,r 

where z is the distance from the internal surface of the glazing, positive in the direction of the external side. 

Equation (28) describes an energy balance imposed at the surfaces of the /th glazing. Equations (29) and (30) 
define the radiosities at the /th glazing, and r ft ,- = 1 - z t - s^j and r bi = 1 - z t - e bi . 

The temperature difference across the /th glazing is given by Equation (31). It is assumed that the solar 
energy is absorbed uniformly through the thickness of the glazing. 



NOTE 



More details regarding Equation (27) to Equation (32) are given in References [5] and [35]. 
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5.3.2 Interaction with the environment 

The effect of boundary conditions imposed by the environment on the window shall be specified. The internal 
and external temperatures, r ftfI+1 and T bQ are: 

^ft,«+1 = ^aUnt (33) 

T b fl = Tai.ex (34) 

The effect of long-wave irradiance at internal and external glazing surfaces is included by setting 

^ft,«+1 = ^gvjnt (35) 

and J b0 = E gyex (36) 

where £ gv int and E„ ex are given by Equations (1 59) and (1 52) in Clause 8, respectively. 
The effect of the convective heat transfer coefficients at the glazing surfaces is included by setting 

"cv,«+1 = "cvjnt \3') 

A cv,1 = ^cv.ex (38) 

5.3.3 Convective heat transfer coefficient — glazing cavities 

5.3.3.1 General 

Convective heat transfer coefficients for the fill gas layers are determined in terms of the dimensionless 
Nusselt number, Nu f : 



h cvi = Nu, 



(39) 



gv,r 



\"QM,i 



where d gv i is the thickness of the fill gas layer (or pane spacing) i, and /L v i is the thermal conductivity of the 
fill gas. Nu ( is calculated using correlations based on experimental measurements of heat transfer across 
inclined air layers. Nu t is a function of the Rayleigh number, Ra t , the cavity aspect ratio, A gv t , and the cavity 
slope, j. 

It should be recognized that deflection of the panes in high aspect ratio cavities can occur. This deflection may 
increase or decrease the average cavity width, d. This deflection can be caused by changes in the cavity 
average temperature, changes in the cavity moisture content, nitrogen absorption by the desiccant or changes 
in the barometric pressure (due to elevation and/or weather changes) from the conditions during assembly. 

NOTE Reference [6] discusses the effects of glass pane deflection and methods to estimate the change in the 

thermal transmittance due to this deflection. 

The Rayleigh number can be expressed as (omitting the "/" and "gv" subscripts for convenience): 

p 2 d 3 g/3c v AT 
Ra = (dimensionless) (40) 
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Treating the fill gas as a perfect gas, the thermal expansion coefficient of the fill gas, /?, is: 



(41 ; 



where T m = fill gas mean temperature in kelvins. 
The aspect ratio of the fill gas cavity i, is: 
H 



A 



gv,i 



■ gv,( 



(42) 



where H is the distance between the top and bottom of the fill gas cavity which is usually the same as the 
height of the window view area. 

Correlation to quantify convective heat transfer across glazing cavities is presented in the 5.3.3.2 to 5.3.3.6. 
Each of these subclauses pertains to one particular value, or range, of tilt angle, y 

This categorization, as a function of y is based on the assumption that the cavity is heated from the internal 
side (i.e., 7ft ; > T bi _^). If the reverse is true (7f U < T bj _^,) it is necessary to seek the appropriate correlation on 
the basis of the complement of the tilt angle, i80°-'y instead of ^and to then substitute 180°- y instead of / 
when the calculation is carried out. 

y= is horizontal glazing, heat flow upwards 

y= 90 is vertical glazing, heat flow upwards 

y= 180 is horizontal glazing, heat flow downwards 

5.3.3.2 Cavities inclined at < y< 60° 



Nit; =1 + 1,44 



1 708 A 



Racos(y) 



J. 



1,6/ 



1708sin''°(1,8r) 



Racos(y) 



Ra< 10 5 and^ gvi >20 



where (x) 



. (x + \x\) 



' RaCOs(y)^ 
5 830 



(43) 



(44) 



NOTE For more details, see Reference [7]. 

5.3.3.3 Cavities inclined at y= 60° 

Nu = (Nu-\,Nu?) 



where Nu 



1- 



0,093 6Ra 
1 + G 



0,314 ^i 



1/7 



(45) 



(46) 



Nil- 



0,104 



0,175 
A Q,i . 



Ra 



0,283 



(47) 
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0,5 



1- 



Ra 



20,6 



0,1 



(48) 



3160^ 
NOTE For more details, see Reference [8]. 

5.3.3.4 Cavities inclined at 60° < y< 90° 

For layers inclined at angles between 60° and 90°, a straight-line interpolation between the results of 
Equations (45) and (49) is used. These equations are valid in the ranges of: 

10 2 <i?a<2x10 7 and5<^ gvr < 100 

NOTE For more details, see Reference [8]. 

5.3.3.5 Vertical cavities 

Nu = (Nui,Nu?) 

Nu^ =0,067 383 8Ra 

Nui = 0,028 154 J Ra°' 4134 

Nu t = 1 + 1, 759 667 8 x 10' u Ra ' 



1/3 



-10 „ 2,298 475 5 



5x10 4 <Ra 

10 4 <^a«:5x10 4 

Ra^^0 4 



Nu 2 =0,242 



Ra 
V A QVJ 



\ 



0,272 



NOTE For more details, see Reference [9]. 

5.3.3.6 Cavities inclined from 90° to 180° 

Gas layers contained in downward facing windows are modelled using: 

Nu = / \ + {Nu^ l -l)sin/ 
Nu v is the Nusselt number for a vertical cavity given by Equation (49). 

NOTE For more details, see Reference [10]. 

5.3.3.7 Fill gas properties 

The density of fill gases in windows is calculated using the perfect gas law. 
PM 



P 



KT n 



P = 1 01 300 Pa and T m = 293 K. 



(49) 
(50) 
(51) 
(52) 

(53) 



(54) 



(55) 



The specific heat capacity at constant pressure, c and the transport properties /u and X are evaluated using 
linear functions of temperature, e.g., the viscosity can be expressed as: 



jU = a + bT n 



(56) 
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Values of a and b coefficients appropriate for calculating c , /u and X for a variety of fill gases are given in 
Annex B. 

5.3.4 Properties of fill gas mixtures 

5.3.4.1 General 

The density, conductivity, viscosity and specific heat capacity of gas mixtures can be calculated as a function 
of the corresponding properties of the individual constituents. 

NOTE For more details, see Reference [11]. 

5.3.4.2 Molecular mass 



!=1 

where m i is the mole fraction of the gas component i, in a mixture of gases. 



(57) 



5.3.4.3 Density 



KT„ 



(58) 



5.3.4.4 Specific heat 



M n 



(59) 



where 

y 

c Pmix = 2X' c P,i 

and the molar specific heat capacity of the ; th gas is 

c p,i ~ c p,i i 

5.3.4.5 Viscosity 



(60) 



(61) 



/"mix _ 2^i 



Mi 



i=1 



7=1 

V J* 1 



m 



111: 



where 



(62) 



>£ 



1/2 



^ + (M i /Mj) (MjIMi) 



J/4 



2^2 \ + (M i lM j \ 



1/2 



(63) 
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5.3.4.6 Thermal conductivity 

^mix = A- mix" 1 "^ mix (°4) 

where 

X' is the monatomic thermal conductivity; 

X" is the additional energy moved by the diffusional transport of internal energy in polyatomic gases. 



•^ mix ~~ / , 



x\ 



;=1 



^ m j 



7=1 



and 



W ■ ■ = 



1 + 



J/2 



(A'M'j) (MilMj) 



1/4 



2v / 2ri + (M,-/M y ) 



1/2 



1 + 2,41 



(m ; — Mj)(M t - 0.142M y ] 

(M i+ Mj) 2 



and 



A" 



mix ~~ / , 
i=1 



X\ 



A m j 



7=1 



where the previous expression for <zS r can also be written as 



1/2, 



+ (X' { /X'j) (m,./m ; ) 



,1/4 



1/2 



^'J : 

2>/2"[i+(jif I -/jif_ / ) 

To find /l mix , use the following steps. 

a) Calculate A',- 

3 , 15 9? 

4 Mf ' 

b) Calculate A",- 

X" t = X t - X' i 

X t is the conductivity of the /th fill gas component (see Annex B) 

c) Use X' t to calculate A' mix 

d) Use X" t to calculate A" mix 



e ) ^mix _ ^ mix +/ ^ mix 



(65) 



(66) 



(67) 



(68) 



(69) 



(70) 



(71) 
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6 Frame effects 

6.1 Area and linear thermal transmittance 

Frame regions of the fenestration system consist of opaque areas that may or may not contain air cavities. 
Frames can be made from a variety of materials, but most common materials are wood, vinyl, aluminium and 
combinations of those (e.g., vinyl-clad wood). 

Frame area thermal transmittance (thermal transmittances) and linear thermal transmittance, X F, shall be 
determined using two-dimensional numerical modelling. This two-dimensional analysis shall provide the rate 
of heat transfer through each unique frame section. See national standards for the required cross sections to 
be considered. Details regarding the required two dimensional numerical analysis are provided in 6.2 to 6.6 
(also see ISO 10077-2). 

6.2 Governing equations for calculating thermal transmittance 

The governing equation shall be developed by imposing an energy balance describing steady-state heat 
transfer by conduction. The governing equation shall be discretized using a conservative formulation (i.e., the 
evaluation of energy flow between two specific nodes or across any given control volume face shall be done in 
a consistent manner throughout the analysis). The frame/edge-glass geometry and the corresponding thermal 
conductivity for each of the various materials, A fe , shall be specified. The numerical solver shall be able to 
generate the two-dimensional heat flow and temperature patterns that satisfy the governing equation. In 
Cartesian co-ordinates, this equation is: 

d 2 T d 2 T q m sc 



(72) 



dx dy *, 

where q\ represents the internal heat generation in watts per square metre. 

The density of heat flow rate, q, is conserved across any surface where two materials meet and is given by: 
AffF dT \ 

q = - l \j-, e ^^-y e y) (73) 

where e x and e are the components of the normal vector to the surface. 

At the boundary, the density of heat flow rate, q, is equal to: 

q = q C v + 1r + ?sc(?sk) ( 74 ) 

where q cv is convective component, and q x is radiative component of the density of heat flow rate, which shall 
be determined in accordance with 8.3 and 8.4 respectively. Quantity q sc (q sk ) is the prescribed density of heat 
flow rate at the boundary (source or sink). 

6.3 Geometric representation and meshing 

6.3.1 Geometric representation 

A two-dimensional representation or model of each frame, sash and edge-glazing assembly shall be made. 
The dimensions of all parts shall be the nominal values as given on the manufacturer's drawings, provided 
that these drawings are an actual representation of the fenestration product. Small radii and minor variations 
in material thicknesses due to manufacturing tolerances or strengthening/attachment requirements may be 
ignored. 

Reinforcing or operating hardware that is essentially continuous, assembly screws or bolts that extend from 
the internal to the external side or bridge a thermal break, including any incompletely de-bridged thermal 
break, shall be included in the model. These thermal-bridging elements may be modelled with three- 
dimensional computational tools when available, otherwise they shall be modelled using the procedure 
outlined below: 
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NOTE For more details, see Reference [34]. 

Calculate the effective conductivity of thermal bridging elements (e.g., bolts, screws, etc.) 

^eff =^br x ^br+(1-^br)^o (75) 

where 

F br = slw 

s is the size of thermal bridging element (e.g., size of a bolt head); 
w is the spacing of thermal bridging elements; 
A br is the conductivity of thermal bridging material; 
X Q is the conductivity of the cross-section without the thermal bridge. 
Use the following criteria to determine if it is necessary to apply the above procedure. 

a) If F br < 1 %, ignore thermal bridge; 

b) If 1 % < F br < 5 %, model using the above method providing that X bx > 1 X ; 

c) If F br > 5 %, always model using the above method. 

Components in the window assembly, which are compressed or deformed from their original shape once 
installed in the window (e.g. weather-stripping), shall be modelled in the compressed or deformed 
configuration. Adjustments to the dimensions of the geometric model are allowed only if they have no 
significant influence on the calculation (see ISO 10211-1). 

Those segments of the actual cross section that are made up of vertical and horizontal surfaces shall be 
represented by similar straight lines that preserve the nominal thickness and relative position of the segment. 
Sloped lines shall be represented by a) similar sloped lines that preserve the nominal thickness and relative 
position of the segment or by b) a series of horizontal and vertical lines, which meet criteria 1 to 4 below. 
Curves shall be represented by a) similar curves that preserve the nominal thickness and relative position of 
the segment or by b) a series of horizontal, vertical, and sloped lines or a series of horizontal and vertical lines 
which meet criteria 1 to 4 below (see Figure 7). 

1 ) The thickness of the representation, d, is equal to the average thickness. 

2) All points on the represented line are within 5 mm of the actual line/curve. The averaged distance (for 
all points) between the represented line and the actual line/curve is less than 2,5 mm. 

3) For conductive materials (materials where the conductivity is 10 times or more than that of any 
surrounding material), the path length shall be maintained to within 5 %. If this condition is not 
possible, then the product of web thickness times conductivity shall be replaced by web thickness 
times conductivity times (cos /+ sin y), where y is the angle of inclination of the sloped web. The 
same result is obtained whether y is measured from the vertical or the horizontal reference. 

4) When sloped materials are represented by a series of rectangles, the contact length between 
adjacent rectangles or polygons, /, is equal to the average actual thickness, d. 

Some windows have nailing flanges that are used to help secure the window in the rough opening. If these 
flanges are intended to be covered up by the exterior cladding (e.g., siding or brick), the portion of the flange 
extending outside the rough opening shall be ignored. 
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In most cases, the internal and external boundaries shall follow the frame profile. In the case where there are 
exterior and interior open frame profile cavities (ventilated cavities and groves) the procedure described in 6.7 
shall be followed. 






.-'' / 




y / 






s / 


r~ 





^ 



z^ 




1 

Key 

1 actual 

2 preferred 

3 acceptable 



Figure 7 — Examples of possible approximations of the actual cross section 



6.3.2 Meshing 

The two-dimensional geometric model shall be divided or meshed into a series of small elements in order to 
provide an accurate representation of the heat flow patterns and temperature distributions. Mesh resolution 
shall be sufficient to ensure that the combined frame/edge thermal transmittance for each cross-section, 
obtained by solving the governing two-dimensional heat transfer equation given in 6.2, shall be within 1 % of 
the combined frame/edge thermal transmittance obtained from an ideal (i.e., infinitely fine) mesh. Acceptable 
meshing schemes include the following: 

a) Successive refinements: The governing heat transfer equation is solved for thermal transmittance for a 
given meshing arrangement. The mesh is made finer either uniformly or in regions of high two- 
dimensional heat flow and a new thermal transmittance determined. An extrapolation is made to the 
thermal transmittance with an infinite number of nodes. The mesh is fine enough when the calculated 
thermal transmittance is within 1 % of the extrapolated thermal transmittance. 

NOTE This requirement is more stringent than that specified in ISO 10211-1, which requires that the number of 
subdivisions be doubled until the change in heat flow through the object is reduced to a prescribed tolerance. The 
more stringent criterion, specified above, is now possible with the increase of computing power. Finite element and 
finite volume methods, with unstructured (non-rectangular) meshes, can also meet this more stringent criterion using 
error estimation methods such as the one given in b) below. 
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b) Energy error norm [13, 14] is applied so that the calculated frame/edge thermal transmittance is within 
1 % of the thermal transmittance determined with an ideal mesh. 

c) Any other approach documented in refereed publications applied so that the calculated frame/edge 
thermal transmittance is within 1 % of the thermal transmittance determined with an ideal mesh. 

6.4 Solid materials 

The thermal conductivity values are usually taken from national standards. Where this is not the case the 
values listed in ISO 10077-2 may be used, but only if they directly match the materials used in the window 
construction. If neither of these sources is used, the thermal conductivity values are to be determined in 
accordance with ISO 8302 (guarded hot plate) or ISO 8301 (heat flow meter) at a mean temperature 
appropriate to national standards. It is assumed that all material thermal conductivity values are constant with 
respect to temperature. 

The surface emissivity values of frame materials are usually taken from national standards. Where this is not 
the case surface emissivity values shall be determined in accordance with ISO 10077-2, but only if they 
directly match the materials used in the window construction. 

6.5 Effective conductivity — Glazing cavities 

Cavities shall be treated as if they contain an opaque solid with an effective conductivity. The effective 
conductivity of a given cavity shall be calculated using the results of the vision area analysis. At cavity i: 



^eff.r 



f r 

gv,( 



b./-1 



(76) 



6.6 Effective conductivity — Unventilated frame cavities 

6.6.1 General 

A frame cavity shall be treated as though it contains an opaque solid, which is assigned an effective 
conductivity. This effective conductivity accounts for both radiative and convective heat transfer and shall be 
determined as follows. 

^eff ={ h cy +h r )xd (77) 

where 

/l eff is the effective conductivity; 

h cv is the convective heat transfer coefficient; 

h r is the radiative heat transfer coefficient {h r = in the case when detailed radiation procedure is 
used); 

d is the thickness or width of the air cavity in the direction of heat flow. 
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The convective heat transfer coefficient, h cv , is calculated from the Nusselt number, Nu, which can be 
determined from various correlations, depending on aspect ratio, orientation and direction of heat flow. 



L„ = Nu- 



(78) 



There are three different cases to be considered, depending on whether the heat flow is upward, downward, 
or horizontal. 



6.6.2 Heat flow downward 

Nu=\0 



(79) 



See Figure 8. 




Key 

1 q = 

2 T ch 

3 T cc 

Figure 8 — Illustration of rectangular frame cavity downward flow direction 
6.6.3 Heat flow upward 

This situation is inherently unstable and will yield a Nusselt number that is dependent on the height-to-width 
aspect ratio, L y IL h , where L y and L h are the largest cavity dimensions in the vertical and horizontal directions, 
see Figure 9. 
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Key 

1 q = 

2 T cc 

3 T cU 



Figure 9 — Illustration of rectangular frame upward flow direction 



a) for — — < 1 convection is restricted by wall friction, and 



AT« = 1,0 



(80) 



b) for 1 < — — < 5 the Nusselt number is calculated according to the method given by 



Nu = ^ + 



Ra 



Jt1 + 2(*2) 



1-lnfc2 



Ra \ 
5 830 



1/3 



-0,95 



I tor 



-1 



(81) 



where 

^1 = 1,40 



Jt2 = 



Ra 



1/3 



450,5 



(82) 
(83) 



W 



X+ X 



(84) 



i?a cr is a critical Rayleigh number, which is found by least squares regression of tabulated values. 
NOTE For more details, see Reference [15]. 
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Ra cr = e 



0,721^+7,46 



(85) 



Ra is the Rayleigh number for the air cavity: 



Ra 



Pa\ L vgP c p,a\( T ch- T cc) 
/"ai^ai 



(86) 



c) for — — > 5 the Nusselt number is: 



AT M = 1 + 1,44 1 



1708 V 

Ra J 



Ra •' 
5 830 



1/3 



(87) 



NOTE For more details, see Reference [7] 

6.6.4 Horizontal heat flow 

See Figure 10. 

2 
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Figure 10 — Illustration of rectangular frame cavity horizontal flow direction 



L 1 
a) for — ^-< — the Nusselt number is: 
L h 2 



Nu = 1- 



2,756x10~ 6 J Ra 2 



i j v 



V L h) 



-0,386 



0,623/Ja 1 ' 5 



f 1 ^ 2/5 



-0,386 



-2,59 



(88) 
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NOTE For more details, see Reference [15]. 

where Ra is Rayleigh number and is defined as: 

.2,3. 



Ra 



Pa\ L hgfi c p,a\{ T ch- T cc) 
/"ai^ai 



(89) 



b) for — - > 5 the following correlation, also the maximum Nu [= (Nu^, Nu 2 , Nu 3 ) max \ is given as: 



Nu 



1 =i 



(0,104^ - 293 



6310 
Ra 



1,36 



1/3 



(90) 



Nu 9 =0,242 



Ra- 



\ 



"v J 



0,273 



(91) 



Nu 3 = 0,060 5Ra 



1/3 



(92) 



NOTE 



1 



For more details, see Reference [9]. 



c) for — < — < 5 the Nusselt number is found using a linear interpolation between the endpoints of a) and 
2 L h 
b) above. 

For jamb frame sections, frame cavities are oriented vertically and therefore the height of the cavity is in the 
direction normal to the plane of the cross section. For these cavities it is assumed that heat flow is always in 
horizontal direction with L y /L u > 5, and so correlations in Equations (90) to (92) in 6.6.4 b) shall be used. 

The temperatures T ch and T cc are not known in advance, so it is necessary to estimate them. From previous 
experience it is recommended to apply T ch = 10 °C and T cc = °C. However, after the simulation is done, it is 
necessary to update these temperatures from the results of the previous run. This procedure shall be 
repeated until values of T ch and T cc from two consecutive runs are within 1 °C. Also, it is important to inspect 
the direction of heat flow after the initial run, because if the direction of the bulk of heat flow is different than 
initially specified, it will need to be corrected for the next run. 

For an unventilated, irregularly shaped frame cavity, the geometry shall be converted into the equivalent of a 
rectangular cavity in accordance with the procedure in ISO 10077-2 (see also Figure 11). For these cavities, 
the following procedure shall be used to determine which surfaces belong to vertical and horizontal surfaces 
of equivalent rectangular cavity (see also Figure 12). 

If the shortest distance between two opposite surfaces is smaller than 5 mm, then the frame cavity shall be 
split at this "throat" region. Also: 

a) any surface whose normal is between 315° and 45° is a left vertical surface; 

b) any surface whose normal is between 45° and 135° is a bottom horizontal surface; 

c) any surface whose normal is between 135° and 225° is a right vertical surface; 

d) any surface whose normal is between 225° and 315° is a top horizontal surface. 
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Figure 11 — Illustration of the treatment of irregularly shaped frame cavities 

b 1 




Key 

a left vertical surface 

b bottom horizonal surface 

c right vertical surface 

d top horizontal surface 

1 direction of the normal to the 
surface 



Figure 12 — Illustration of how to select surface orientation for frame cavities; dashed lines indicate 
direction of the normal to surface with cut off angles at 45°, 135°, 225° and 315° 
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Temperatures of equivalent vertical and horizontal surfaces shall be calculated as the mean of the surface 
temperatures according to the classification shown above. The direction of heat flow shall be determined from 
the temperature difference between vertical and horizontal surfaces of the equivalent cavity. The following rule 
shall be used (see also Figure 13): 

a) heat flow is horizontal if the absolute value of the temperature difference between vertical cavity surfaces 
is larger than between horizontal cavity surfaces; 

b) heat flow is vertical, heat flow upwards if the absolute temperature difference between horizontal cavity 
surfaces is larger than between vertical cavity surfaces and the temperature difference between the top 
horizontal cavity surface and bottom horizontal cavity surface is negative; 

c) heat flow is vertical, heat flow downwards if the absolute temperature difference between horizontal cavity 
surfaces is larger than between vertical cavity surfaces and the temperature difference between the top 
horizontal cavity surface and bottom horizontal cavity surface is positive. 




Key 

1 temperature of bottom surface, T bo 

2 temperature of left vertical surface, 7j v 

3 temperature of top surface, T ip 

4 temperature of right vertical surface, T u 

a ) Krt-^lfp Ftp-^bo 



b) |r„ - 7| f | < r tp - r b0 and r tp < r b0 



c) |r rt -r| f |<|r tp -r b0 | and r tp >r, 



bo 



heat flow is horizontal; 



heat flow is vertical, heat flow up; 



heat flow is vertical, heat flow down. 



Figure 13 — Illustration of how to select heat flow direction 
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6.6.5 Radiant heat flow 

The radiative heat transfer coefficient h r shall be calculated using: 



h r 



1 1 

— + 



4 a T~ 



1 



1 + 



\ L M J 



Lu 



(93) 



where T av = -^ — 



(94) 



The above notation assumes radiant heat flow to be in the horizontal direction. If the heat flow direction is 
vertical, then the inverse of the ratio L h /L v shall be used (i.e., LJL h ). 

NOTE For more details, see Reference [16]. 



6.7 Ventilated air cavities and grooves 

6.7.1 Slightly ventilated cavities and grooves with small cross section 

Exposed grooves with small cross sections (see Figure 14) or cavities connected to the external or internal 
environments by a slit greater than 2 mm but not exceeding 10 mm are to be considered as slightly ventilated 
air cavities. The equivalent conductivity is twice that of an unventilated air cavity of the same size in 
accordance with 6.6. For cases where the slit is less than or equal to 2 mm, treat the cavity as completely 
enclosed in accordance with 6.6. 

Dimensions in millimetres 
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Figure 14 — Examples for slightly ventilated cavities and grooves with small cross section 
6.7.2 Well-ventilated cavities and grooves with large cross section 



In cases not covered by 6.6 and 6.7.1, in particular when the width b of a groove or of a slit connecting a 
cavity to the environment exceeds 10 mm, it is assumed that the whole surface is exposed to the 
environment. Therefore, the surface heat transfer coefficients, h m \ and h ou{ , calculated in accordance with 
Clause 8, shall be used at the developed internal and external surfaces, respectively. 
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In the case of a large cavity connected by a single slit and a developed cross sectional length exceeding the 
width of the slit by a factor of 5 (see Figure 15), the detailed radiation model (see 8.4.2) may be appropriate 
for determining the radiative portion of the surface heat transfer coefficient. 
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Figure 15 — Examples for well-ventilated cavities and grooves 
7 Shading devices 

7.1 Definitions 
7.1.1 Introduction 

Clause 7 provides the necessary equations for the effects of shading devices on the thermal and optical 
properties of a window system. 

The scope is restricted to those kinds of shading devices which are or which may, by proper approximation, 
be treated as a layer parallel to the pane(s) of the window. 

The introduction of shading devices in the model of the window system leads to modifications of the main 
thermal and optical equations. In order not to complicate the presentation of equations given in Clause 5, the 
necessary changes to those equations are not integrated in Clause 5 itself, but given as amendments in 
Clause 7. 

NOTE Information on calculation procedures and measurement techniques on shading devices can be found in 

references [17 - 25]. In general, these references concern ongoing work. 

The contents of this clause are based on the most up-to-date procedures, with simplifying approximations 
where needed due to practical limitations with respect to modelling and computational efforts and availability 
of product data. 

Shading devices can be divided into two basic types: 

— layer type of shadings, such as screens, curtains and Venetian blinds which are located parallel to the 
pane(s), with intimate thermal-optical contact; 

— extra-fenestrial type of shadings, such as awnings and overhangs, which are located less close to the 
panes. These types of shading may be regarded as part of the window's environment, because of the 
limited thermal interaction. They have mainly an effect on the temperature and radiation conditions 
outside the window. In specific cases, however, the conditions and properties of the window itself may 
also influence the condition of this "environment" (e.g. reflection of solar radiation, hot air pockets). These 
types of shading are outside the scope of this International Standard. 
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7.1 .2 Principle of the calculations 

The thermal-optical interaction of a layer type of shading device is, to a greater extent, similar to the panes 
and films. In this regard, the layer type of shading device may be defined in the model as a layer between two 
gaps. This thus-defined layer exchanges heat with the other components and/or the environment by 
conduction and convection and by thermal radiation. It also absorbs, reflects and transmits solar radiation. 
However, due to its porous structure (open weave, slats), the shading device is not only partially transmittant 
for solar radiation, but also for thermal (long wave) radiation. It shares this characteristic with some suspended 
thin films. This phenomenon is already covered in the equations by introducing in the equations transmittance 
for thermal radiation. 

The shading device is usually also permeable for air, either due to its porous structure or due to openings at 
its perimeter. Air may cross the shading device and thus move from one gap to the other or from the 
environment into the gap behind the shading device and vice versa. This phenomenon has not been 
previously covered by the equations in the previous clauses and will therefore be introduced in this clause. 

Because the layer type of shading device is modelled as a one-dimensional layer similar to a pane or film, the 
two- or three-dimensional characteristics shall be translated into one-dimensional numbers. This is in particular 
the case for the optical properties, e.g., the optical properties of a shading device are a function of the geometry 
of the device and the position in the assembly. To consider a slat type of shading device such as a Venetian 
blind, information on the optical properties of the slat material, together with the geometry of the slats and their 
positions, is used to determine the overall transmittance, reflectance and absorption of the layer. 

7.2 Optical properties 

7.2.1 General 

A particular characteristic of a shading device compared to "normal" glazings or films is that the incident solar 
radiation may change direction while being transmitted or reflected at the layer. 

For the evaluation of thermal effects, the following approximation is considered to be sufficiently accurate. 

Beam radiation transmitted or reflected by the solar shading device is considered to be split into two parts: 

— an undisturbed part (specular transmission and reflection); 

— a disturbed part. 

The disturbed part is approximated as anisotropic diffuse (Lambertian). 

Diffuse radiation transmitted or reflected by the solar shading device is assumed to remain diffuse. 

An exact description of the way solar radiation travels through the system would require a full three- 
dimensional calculation using the full matrix of the transmission, absorption and forward and backward 
reflection for each angle of incidence at each component. For the evaluation of the spatial distribution of 
daylighting, this would be the necessary way to proceed. 

Consequently the following solar properties of the solar shading device are required or transmittance, for 
beam radiation, for each angle of incidence: 

r dir dir^w/) direct to direct transmittance; 

r dir,difO*vitf) direct to diffuse transmittance. 

For diffuse radiation: 

r dif diK^w/) diffuse to diffuse transmittance. 
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Similarly for the reflectance, the following properties are required (for beam radiation, for each angle of 
incidence): 

r dir dir(^w/) direct to direct reflectance; 

r dir dif(^w/) direct to diffuse reflectance. 

For diffuse radiation, the following properties are required: 

r dif diK^w") diffuse to diffuse reflectance, and for the absorption: 

a 4\AKij) = H " r dir,dir( /l w/) - r dir,dir( /l w/) ~ ^dir.diKA^) ~ ^ir.dif^w/)] ( 95 ) 

«dif(^w/) = [ 1 " ^dif,dif(-V) " r dif,dif(-V')l (96) 

7.2.2 Resulting amendments to the equations in 5.2 and 5.3 

For a window system incorporating a layer type of shading device, the optical equations given in Clause 5 
remain the same with the following extensions. 

a) Each spectral flow rate equation in 5.2 and 5.3 shall be split into three: a "dir,dir", "dir,dif" and "dif,dif" heat 
flow, with corresponding transmittance rand reflectance r. In the sum of the spectral heat flows: "dir,dir", 
"dir.dif" and "dif.dif" parts shall be summed. 

b) The transmittance shall be split, similar to the reflectance, into a forward and a backward value. 

c) The sum of r dird j r and r dirdif is equal to the direct to hemispherical transmittance r dirh ; similarly for the 
reflectance. 

For slat type shadings, equations to calculate these properties are given in 7.3, on the basis of optical 
properties and the geometry of the slats. 

NOTE There is no existing International Standard for the measurement of these optical properties. Until such a 

testing standard is available, the calculation method of this clause shall be considered as provisional and is provided for 
information purposes only, except for slat types of shading devices for which the following clause provides a calculation 
method. 

Once a beam transmitting through or reflecting at a solar shading device is split into a direct and a diffuse part, 
the diffuse part continues its route through the system. This implies that even for normal incidence solar 
radiation, for all other panes, films and shading layers in the window, the r difdif and r difdif values are required; 
consequently, the values for normal incidence provide insufficient information. 

Due to the redirection of the radiation, the forward transmittance is not necessarily equal to the backward 
transmittance, as illustrated in Figure 16. 
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Figure 1 6 — Illustration of different values for forward and backward solar transmittance 
(slats with different colour at both surfaces) 



7.3 Slat type of shading 

7.3.1 General 

For a shading device consisting of parallel slats, the optical properties can be determined as function of slat 
properties, geometry and position (see Figure 17). The air permeability can also be determined as function of 
slat geometry and position. 
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Figure 17 — Slat geometry 
7.3.2 Optical properties 

7.3.2.1 General 

This subclause gives the procedure for calculating the solar optical properties of a slat type of shading device 
provided: 

— that the slats are non-specular reflecting; 

— that any effects of the window edges may be ignored. 

The procedure is to consider two adjacent slats and to subdivide the slats into five equal parts (see Figure 1 8). 
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Figure 18 — Discretization used in the model 



Every slat is divided into five elements (the improvement of considering more elements is negligible). Notice 
that different properties can be assigned to every element, in particular to every side of the slat. The process 
described below has to be solved for every wavelength band required by the properties of the elements or by 
the rest of the transparent system where the shading device is installed. 

Due to the assumption of non-specular reflection, a slight curving of the slats may be ignored. 

7.3.2.2 Equations 

These equations have a more general application, if the allocation of the layer numbers is generalized. 
For each layer f,z and b,i, with i from to n (here: n = 6) and for each spectral interval A^,-, (^->^ + A/l^,): 

E U = X[( r W + T b,k ) E f,k F u^,i + { r b,k + *u ) E b,k F b,k^,i] ( 97 ) 



E b,i = X[( r W + r U ) E b,k F b,k^b,i + {n,k + *b,k ) E f,k F f,k^b, 



(98) 



where 



E k is the irradiance on surface k; 



^p^.q is the view or shape factor from surface p to surface q (e.g. b,k to b,0- 
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Also 

^f,ex = "'exlAw/v 

^b.int = J \r&(Ktj) = ° 

where 

J eK is the radiosity from the external environment (incident solar radiation); 

J int is the radiosity from the internal environment (room reflection). 



(99) 
(100) 



7.3.2.3 Diffuse-diffuse transmission and reflection 

Due to the assumption of non-specular reflection, the values for the view factors ^p_>q can be calculated by 
conventional view factor calculation methods for diffuse radiation exchange. 

NOTE For calculation methods on view factors, see [26]. 

For diffuse incident radiation, the view factor between external environment and the other layers is also 
determined by the view factors for diffuse radiation exchange. 

After solving the set of equations, the diffuse/diffuse transmission coefficient is found as the radiation reaching 
the internal environment E fn (n = 6), divided by the incident solar radiation, J ex : 



r dif,dif (^w j ) - E f,n (^w j ) I J ex {^w j ) 
Similarly, for the diffuse/diffuse reflection coefficient: 

r dif,dif (^w j J = E b,ex (^w j y J ex (^w j ) 



(101) 



(102) 



7.3.2.4 Direct-direct transmission and reflection 

By straightforward geometric calculation from the angle and aspect ratio of the slats (see Figure 19), the beam 
radiation which passes the slats without touching can be calculated for given angle of incidence y. 





Figure 19 — Direct-direct transmission 



©ESA 



37 



ESISO15099:2012(E) 

This part of the transmission is wavelength-independent. 
This is the direct/direct transmission: E dir d - ir (/ w ) 

Consequently, the direct-direct transmittance for incidence angle /amounts to: 

Z"dir,dir {/) = ^dir.dir (^wy.rj/^ex ( A vij>r) (103) 

for any wavelength A^-. 
There is no reflected radiation to the outside without reflecting on one or more parts of the shading device, so: 

''dir.dir (r) = (104) 

7.3.2.5 Direct-diffuse transmission and reflection 

Firstly, calculate, for the given angle of incidence y, which parts of the shading device k are directly irradiated 
by J f ex (see Figure 20). 




Key 

1 outdoor 

2 indoor 

Figure 20 — Directly irradiated parts of the shading device 

The view factors between the incident radiation J ex and those directly irradiated parts k are: 



^f,ex->f,* - 1 



^f,ex->b,jfc ~ 1 
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The view factor between the internal and external environment is zero, to exclude the direct-direct 
transmittance: 

F f,ex->b,« = ° and F b,ex->f,H = ° 
After solving the set of equations we find the direct-diffuse transmission and reflection coefficients: 

z"dir,dif [Kj/] = E U {A Wj ,r)lJex (^wy./J ( 105 ) 

r dir,dif \KjiY) = E b,n {Awj>r)lJ e x ( A w 7 .r) (106) 

7.3.2.6 Absorption 

That part which is neither transmitted, nor reflected, is the part which is absorbed in the slats. Per wavelength 
band: 

(107) 
(108) 



"dir (^w j)~\ 1 ~ r dir,dir (^w j ) ~ r dir,dir (^w j j ~ r dir,dif (^w j j ~ r dir,dif (^w j j 
«dif (^w j j = | 1 ~ ^dif.dif (^w j ) ~ Pdif.dif (^w j j 



7.3.2.7 Thermal transmittance and reflectance 

The blinds are also semi-transparent for infrared (thermal) radiation. In order to obtain the IR transmittance 
and reflectance of the shading device for given (IR) slat properties, the same model is used as for the 
calculation of diffuse-diffuse transmission and reflection of solar radiation, replacing the slat's solar optical 
properties by its thermal radiation properties. 

The normal emissivity of the surfaces can be measured in accordance with EN 12898. There is no existing 
standard for the measurement of the hemispherical emissivity of opaque materials. Usually an emissometer is 
used for this purpose. 

Selected examples of calculated optical properties of a slat type of shading device as function of the slat 
properties and geometries are given in Annex C. 

7.4 Ventilation 

7.4.1 General 

For a ventilated cavity, the set of equations given in 5.3 shall be extended in the way described in 7.4.2. 

7.4.2 Main heat balance equations 

7.4.2.1 Principle 

Air spaces may be connected to the exterior or interior environment or to other spaces. For a ventilated gap, 
the heat balance in the gap requires an extra term, the amount of heat supplied to or extracted from the gap 
air. This implies that it is no longer sufficient to describe, as in 5.3, the conductive/convective heat exchange in 
a gap as the heat transfer from one surface to the other. It is necessary to make a split between the 
conductive/convective heat transfer from one surface to the air and from the air to the other surface as 
illustrated in Figure 6. In the heat balance equations for the gap, the heat extracted from or supplied to the gap 
by ventilation is added to this air gap node. 
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The mean temperature of the air in the gap is given by equations for the heat exchange between the air 
flowing through the gap and the adjacent surfaces. 

NOTE There is no existing International Standard for measurement of these properties. Until such a testing standard 

is available, the calculation method in 7.4.2 must be considered as provisional and is provided for information purposes 
only. 

7.4.2.2 Equations — non-vented gap 

For the non-vented case (5.3) the heat exchange by conduction/convection across a gap from one layer to the 
adjacent layer (pane, film or shading device) as given in 5.3.1: q cv j =h cvi (T^ i -T bi+ ^), is split into two 
parts (see Figure 21 ), with the mean temperature of the air in the gap as a variable. 
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Figure 21 — Split convective heat transfer across non-vented gap 



9cv,f,/ ~~ ^"cdv.i \1 f,i * gap,; J ~~ <7cv,b,/+1 ~~ ^"cdv.r \J gap,r ■* b,;'+1 j 



(109) 



where 

q CV fj is the convective heat transfer from the one surface to the gap, in watts per square metre; 

h^yi is the surface-to-surface heat transfer coefficient by conduction/convection for non-vented 
cavities, given by the equations in 5.3, in watts per square metre kelvin; 



hi 



gap,/ 



is the temperature of the surface of layer (pane, film or shading) i, facing cavity i, see 5.3, in 
kelvins; 



is the equivalent mean temperature of the air in cavity i, given in Equation (119) below, in kelvins; 

lev b 1+^ ' s * ne convective heat transfer from the gap to the other surface, in watts per square metre; 

T b i+1 is the temperature of the surface of layer (pane, film or shading) i + 1 , facing the cavity i, see 5.3, 
in kelvins. 



7.4.2.3 Ventilated gap 

In a ventilated gap, due to the air movement, the convective heat exchange coefficient is increased 
(see Figure 22) 



40 



©ESA 



ES ISO 15099:2012(E) 



\, 



?4j 



gap, /, out 



' f, / + 1 



gap,' 



/7 

CV, I 



:h 



?v./ 



7" 




gap, /, inl 



Key 

1 pane or shading 

Figure 22 — Model of mean air- and outlet temperature and main dimensions 

This increased coefficient is written as h cvi : 

1c\/,b,i = "cdv.z'Ub,; ~~ lgap,i) < anc " <7cv,f,/+1 = "cv./V^gap,; — -'f.z'+l) 
with /z cv ,- given by the equation: 



(110) 



h CV,i = 2fl CV,i + 4V i 



(111) 



where 



q csl b ,- is the convective heat transfer from the one surface to the gap, in watts per square metre; 

h cvj is the surface-to-air heat transfer coefficient by conduction/convection for vented cavities, given 
by Equation (111), in watts per square metre kelvin; 

q cv f (+1 is the convective heat transfer from the gap to the other surface, in watts per square metre; 

h cAvi is the surface-to-surface heat transfer coefficient by conduction/convection for non-vented 
cavities, given by the equations in 5.3, in watts per square metre kelvin; 

V t is the mean air velocity in the gap, see 7.4.4, in metres per second; 

and with (same as for the non-vented case): 

T„ api is the equivalent mean temperature of the air in cavity i, given by Equation (119) in 7.4.3, in 
kelvins; 

T b i is the temperature of the surface of layer (pane, film or shading) i, facing cavity i, see 5.3, in 
kelvins; 
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r f i+1 is the temperature of the surface of layer (pane, film or shading) i + 1 , facing cavity i, see 5.3, in 
kelvins. 

For zero velocity, the equations for the ventilated gap reduce to the equations for the non-vented case. 

Due to the ventilation, an extra term is added to the heat balance equations of the gap given in 5.3. The extra 
term is: 

9vl,/ = Pi xc p x £>vl,i (^gap./jnl ~ T gap,i,out) I { H i L i) C 112 ) 

Equations (114) to (117) are formulated in such a manner to satisfy the following energy balance equation: 

tfvl.i = 9cv,b,r +9cv,f,;+1 C 113 ) 

where 

g v i ; is the heat transfer to the gap by ventilation, in watts per square metre; 

p t is the density of the air in cavity k at temperature 7g ap ■ , in kilograms per cubic metre; 

c is the specific heat capacity of air, in joules per kilogram kelvin (i.e.: 1008); 

<p^\ i is the air flow rate in cavity /, see 7.4.4, in cubic metres per second; 

r gap t in | is the temperature at the inlet of the gap, in kelvins. 

The heat transfer is normalized to 1 m 2 of the aperture area. 

The value of T'gap / ini depends on where the air comes from: either the internal or external air temperature or 
the outlet temperature T gap k out °f * ne 9 a P ^ w '* n which gap i exchanges air; 

T„ ap i out is the temperature at the outlet of the gap, see Equation (1 14), in kelvins; 
L { is the length of cavity i, see Figure 22, in metres; 

H i is the height of cavity i, see Figure 22, in metres. 

7.4.2.4 Heat transfer to internal environment 

The heat transfer to the internal environment shall be extended in a similar way with a term q vU for the heat 
transfer by ventilation by air coming from cavity k. 

Following the convention from Clause 5, with i = n for the internal environment, for all cavities k with air flow to 
the internal environment, n\ 

<lv\,n =^ l Pi- c p-<Pv\,i( T gapJ,ou{- T a\,n) l ( H i xL i) ( 114 ) 

i 

The heat transfer is normalized to 1 m 2 of the aperture area, 
where 

p t is the density of the air in cavity i at temperature r gap r , in kilograms per cubic metre; 

c is the specific heat capacity of air, in joules per kilogram kelvin (i.e.: 1008); 
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<Nu 



is the air flow rate in cavity i, see 7.4.4, in cubic metres per second; 



^gap/out ' s * ne temperature of the air at the outlet of the gap from where the air originates, 
see Equation (116), in kelvins; 

r ai is the indoor air temperature, in kelvins; 



h 



H; 



is the length of cavity i, see Figure 22, in metres; 
is the height of cavity i, see Figure 22, in metres. 



7.4.3 Temperatures in the cavity 

Assuming the mean velocity of the air in the space is known (see 7.4.4), the temperature profile and the heat 
flow can be calculated by a simple model. Due to the air flow through the space, the air temperature in the 
space varies with height (see Figure 23). 
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Figure 23 — Air flow in the gap of a window system 
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The temperature profile depends on the air velocity in the space and the heat transfer coefficient to both 
layers. The air temperature profile in space i is given by: 

-'gap,; \ x j = ^av,/ ~\*av,i ~*gap,i,\n) e \''v) 

where 

Tgap t (x) is the temperature of the air in gap i at distance x from the inlet, in kelvins; 

H Qi is the characteristic height (temperature penetration length), see Equation (117), in metres; 

r gap i in is the temperature of the incoming air in gap i, in kelvins; 

T av { is the average temperature of the surfaces of layers i and i + 1 , given by equation: 

^.fik^tw) (116) 

where 

T bj is the temperature of the surface of layer (pane, film or shading) i, facing cavity i, see 5.3, in 
kelvins; 

7f !+1 is the temperature of the surface of layer (pane, film or shading) i+ 1, facing cavity i, see 5.3, in 
kelvins. 

The characteristic height of the temperature profile is defined by: 

^0,1=^—^ xV t (117) 

where 

H Q ■ is the characteristic height (temperature penetration length), in metres; 

p i is the density of the air at temperature T^- , in kilograms per cubic metre; 

c is the specific heat capacity, in joules per kilogram kelvin (i.e.: 1008); 

bj is the width of the cavity i, in metres; 

/i cv i i is the heat transfer coefficient for ventilated cavities, see Equation (111), in watts per square metre 
kelvin. 

V i is the mean velocity of the air flow in the cavity i, see 7.4.4, in metres per second; 
The leaving air temperature is given by: 

' gap,r,out = * av,i _ y- av,; ~* gap,/,inl ) e ' (11 o) 

where 

r gap i out is the temperature of the air at the outlet of gap i, in kelvins; 
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T avj is the average temperature of the surfaces of layers i and i- 1, given by Equation (116), in 

kelvins; 

T„ ap ■ in | is the temperature of the incoming air in the cavity i, in kelvins; 

H Q i is the characteristic height (temperature penetration length), given by Equation (117), in metres; 

H t is the height of space i, in metres. 

The thermal equivalent (average) temperature of the air in the space i is defined by: 

-'gap,; = ~rr J -'gap,/ \ x ) x " x = *av,i ~ [* gap.i.out ~ * gap.ijnl J (' '") 

' ' 

where 

r ga p i is the equivalent mean temperature of the air in the cavity i, in kelvins; 

T avj is the average temperature of the surfaces of layers i and i + 1, given by Equation (116), in 

kelvins; 

H Qi is the characteristic height (temperature penetration length), given by Equation (117), in 

metres; 

H { is the height of space i, in metres. 

r gap i out is the temperature of the air at the outlet of the gap i, in kelvins; 

r gap j in | is the temperature of the incoming air in gap i, in kelvins; 

7.4.4 Air flow and velocity 

7.4.4.1 Forced ventilation 

If the air flow within the air layer has a known value (for example due to mechanical ventilation), the equations 
given in 7.4.2 and 7.4.3 shall be applied as such, with the air velocity (m/s) given by: 

V t =-p*j- (120) 

where 

V i is the mean velocity of the air flow in cavity i, in metres per second; 

<p v i { airflow rate in cavity /, in cubic metres per second; 

bj the width of cavity i, in metres; 

L t is the length of the cavity i, see Figure 22, in metres. 

NOTE gx,\ d is the air flow rate for the whole area. 
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7.4.4.2 Wind-induced ventilation 

For external shading devices, the thermally-induced ventilation is mixed with wind-induced ventilation. 

The heat exchange by ventilation between the shading and the next layer (pane) can be described on the 
basis of an appropriate value for the air flow or velocity. An appropriate value is to be determined on the basis 
of experiments or calculations [computational fluid dynamic (CFD) modelling]. 

For conservative design calculations, one may treat the cavity flow as forced convection (see 7.4.4.1 ), with the 
value for the air velocity V; set to extreme low and extreme high values, thereby giving two values for the total 
solar energy transmittance. 

7.4.4.3 Thermally-driven ventilation 

The velocity of the air in the space, caused by the stack effect depends on the driving pressure difference and 
the resistance to the air flow of the openings and the space itself (see Figure 24). 

gap, ;',out gap, /r,in 
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Figure 24 — Schematic presentation of the stack effect 

NOTE The height of the neutral zone x Q depends on the flow resistances of the inlet and outlet openings 

The air velocity is known by solving the set of equations given in 7.4.4. 

The pressure difference results from a temperature difference between the space/ and the connected space 
k, which is the exterior air, the interior air or another space. The temperature profile in the spaces is 
represented by the thermal equivalent temperature [Equation (119)]. The driving pressure difference Ap T may 
be written approximately as: 



^Pj,i,k = PQxTQXgxHjX \cosn 



[\ T gap,i ^gap.^lj 
[ T gap,i x T gap,k ) 



(121) 
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where 

AP T; -£ is the driving pressure difference between space i and space k, in pascals; 

p Q is the density of air at temperature T Q , in kilograms per cubic metre; 

r is reference temperature, e.g., T Q = 283 K; 

g is accelaration due to gravity = 9,81 metres per second per second; 

y i is the tilt angle of space i in degrees from vertical; 

H { is the height of space i (same as space k), in metres; 

T„ ap ■ is the equivalent (mean) temperature of the air in the space i, see Equation (119), in kelvins; 

r gap k is the equivalent temperature of the connected space, which may be another gap k or the internal 
or external environment, in kelvins. 

The air flow in the space is described as a pipe flow. Therefore, the following effects have to be taken into 
account. 

Acceleration of the air to the velocity V (Bernouilli's equation): 

AP Bi =^- Vi 2 (122) 

Steady laminar flow (Hagen-Poiseuille law): 

AP HP , ; =12-// r ^- (123) 

Pressure loss in the inlet and outlet openings: 

^Pz,i=^-Vi 2 {Z- m +Z ouU ) (124) 

where 

AP Bi is the Bernouilli pressure loss in space i, in pascals; 

p t is the density of air at temperature r gap , , in kilograms per cubic metre; 

V t is the mean velocity of air flow in cavity i, to be solved with Equation (120), in metres per second; 
(same for k)\ 

AP HP ( - is the Hagen-Poiseuille pressure loss in space i, in pascals; 

ju t is the dynamic viscosity of air at temperature T gap ■ , in pascal seconds; 

H i is the height of space i, in metres; 

bj is the width of cavity i, in metres; 

Z { is the pressure loss factors Z of cavity i, according to Equations (126) and (127). 
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The same equations apply to space k, where V k = V i -b i lb k . 



If the space k is the exterior or interior, V k = is assumed, in which case the pressure loss terms AP Bk and 
AP HP k are zero as well as AP zk , where 



AP 



Z,i,k 



is the pressure loss Z between space i and k, in pascals. 



The total pressure loss shall be equal to the driving pressure difference and this results in the velocities V i and 
V k by solving the equation: 



A^T.a = ^BJ + AP H P,« + AP Z(I - + AP Zik + AP B ,k + ^HP,k 
where 



(125) 



AP Jik is the driving pressure difference between space i and space k, according to Equation (121), in 
pascals; 

AP Bj is the Bernouilli pressure loss in space i, according to Equation (122), in pascals; 

AP HP ; is the Hagen-Poiseuille pressure loss in space i, according to Equation (123), in pascals; 

AP Zi is the pressure loss Zat the inlet and outlet of space i, according to Equation (124), in pascals; 

AP Z k the same as AP Z i but for space k; 

AP B k is the Bernouilli pressure loss in space k, according to Equation (122), in pascals; 

AP HP k is the Hagen-Poiseuille pressure loss in space k, according to Equation (123), in pascals. 

The pressure loss factor, Z, for openings may be estimated from the ratio of the equivalent area of an opening 
A eq to the cross section of the space ^ s (see Figure 25). 




Key 

1 front view 

2 side view 



Figure 25 — Openings in a ventilated gap 
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z 



inl 



0.6x^eq,inl,r , 



(126) 



and 

\2 



■^out 



^S,l ,, 



Q 6x 4jq,out,i 



(127) 



where 

^ s ( . is the cross section of space i; A s ,- = 6,-L,-, in square metres; 

& ; - is the width of cavity i, in metres; 

L { is the length of cavity i, in metres; 

^eqjni,/ ' s tne equivalent inlet opening area of cavity i, according to Equation (128) or (129), in square 
metres; 

^eqoutr ' s * ne equivalent outlet opening area of the cavity i, according to Equation (128) or (129), in 
square metres. 

If the temperature T gap i {T gap ^ of the cavity i (k) is higher than the temperature of the connected space k (i): 

1 -^tp 1 Ahn 

^eqjnl = Ao + 7J "j —7- Mlf + ^rt + ^ho) ^eq.out = ^tp + 77 ~ ^~ (^lf + ^rt + ^ho) ( 128 ) 



2 ^bo+^tp 2 ^bo+^tp 



Otherwise: 



1 ^tp 1 ^hn 

^eq.out = ^bo + 77 —— (4f + ^rt+Ao) ^eqjnl = ^tp + 77 "j 7" Mlf + ^rt + Ao) ( 12 9) 

2 ^bo+^tp 2 ^bo+^tp 

where 

^ s is the cross section of the space, in square metres; 

A ho is the area of the bottom opening, in square metres; 

A ip is the area of the top opening, in square metres; 

A n is the area of the left side opening, in square metres; 

A A is the area of the right side opening, in square metres; 

A ho is the total area of the holes in the surface (homogeneously distributed holes), in square metres. 
It is assumed that the side openings are distributed evenly from top to bottom. 
All these areas are total flow areas for the window (i.e., not normalized). 
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7.4.5 Gas-filled cavity with air circulation 

In those cases of a closed cavity containing a gas mix and other component of the fenestration, e.g., an 
incorporated blind, the gas mix may flow from one side of the component (blind) to the other. In that case the 
equations given above remain valid, if "air" is replaced by "gas-mix", with the corresponding gas-mix 
properties. 

7.4.6 Air permeability of slat types of shading devices 

The air permeability of slat types of shading devices can be described using an appropriate value for the 
equivalent air permeability of the surface, A ho . An appropriate value shall be determined on the basis of 
experiments or calculations [computational fluid dynamic (CFD) modelling]. 

For conservative design calculations, the value for the equivalent air permeability of the surface, A ho , can be 
set to extreme low and extreme high values, thereby giving two extremes for the total solar energy 
transmittance. 

7.5 Total solar energy transmittance and thermal transmittance 

The thermal transmittance, U, and the total solar energy transmittance, r s , shall be calculated using 
Equations (2) and (14) respectively, applying the amendments to the equations in 5.2 as given in 7.2. 



8 Boundary conditions 

8.1 General 

The various thermal properties can be determined using a standard calculation method but each will also be 
affected, to some extent, by the boundary conditions to which the product is exposed, i.e., the environment. 

The boundary conditions consist of: 

a) internal and external air temperatures, T mt and T ex , respectively; 

b) internal and external surface convective heat transfer coefficients, h cv int and h cv ex , respectively; 

c) solar spectral irradiance distribution, E(A V/ ), and a function describing the photopic response of the eye, 
R{A^). Both E(A W ) and R{A^) consist of a set of function values listed for a set of discrete wavelength 
values. Function values at intermediate wavelengths can be found by linear interpolation; 

d) the longwave irradiance on the external and internal glazing surfaces, G g ex and G„ int , respectively, as 
well as the longwave irradiance at the external and internal frame surfaces, G fex and G f int , respectively. It 
is assumed that external longwave irradiance depends on the clearness of the sky factor, f cW . 

8.2 Reference boundary conditions 

8.2.1 General 

Unless a specific set of boundary conditions is of interest (e.g., to match test conditions, actual conditions or to 
satisfy a national standard), the following standard boundary conditions shall be used. In each case the 
following spectra shall be used. 

E S (A^) = ISO 9845-1 (hemispherical solar spectral irradiance tabulated at N s values of X^); 

E y {XJ) = ISO/CIE 1 0526 (colorimetric illuminance tabulated at N sv values of /LJ; 

R{A^) = ISO/CIE 1 0527, (photopic response for the 2° observer tabulated at N v values of 2 W ). 
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8.2.2 Winter conditions 



r int = 20°c 



r ex = o°c 



A CVJnt = 3,6W/(m2.K) 

Vex = 20W/(m2.K) 



T = T 
'r,m *ex 

I s = 300 W/m 2 

8.2.3 Summer conditions 

r int = 25°C 

r ex = 30°c 

/* cvint = 2,5W/(m 2 -K) 
Vex = 8W/(m2.K) 
T = T 
/ s = 500 W/m 2 

8.3 Convective heat transfer 

8.3.1 General 

Convection heat transfer is energy transfer between a surface and a moving fluid. Heat is transferred by 
natural convection (i.e., convection driven by temperature gradient) when the air velocity is sufficiently small 
(i.e., less than 0,3 m/s). On the other hand, heat is transferred by forced and mixed convection for velocities 
above 0,3 m/s. Accurate determination of this convective heat transfer on both internal and external boundary 
surfaces is extremely difficult and can only be done by careful measurements and computer simulation. For 
these reasons, surface heat transfer coefficient correlations had been developed and are given in 8.3.2 and 
8.3.3. 

8.3.2 Convective heat transfer coefficient — internal side 

8.3.2.1 General 

The convective heat transfer on the internal side primarily occurs by natural convection, and rarely by mixed 
and forced convection. Standard boundary conditions assume natural convection on the internal side. The 
density of convective heat flow on the internal boundary is defined as: 



I cv.in 



Vin^s.in-Tin) (130) 



where T s in is the temperature of any internal fenestration surface [i.e., T b n (temperature of the internal glazing 
surface) or the temperature of the internal frame surface]. The convective heat transfer coefficient, h cvin is 
determined from heat transfer correlations given in 8.3.2.2. 
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8.3.2.2 Heat transfer by natural convection 

The natural convection heat transfer coefficient for the internal side, h 
Nusselt number, Nu. 



cv.int' 



is determined in terms of the 



Vint = Nu 



H 



where X is the thermal conductivity of air. 



jt 



(131) 



Nu is calculated as a function of the corresponding Rayleigh number based on the height, H, of the 
fenestration system, Ra H . 



Ra 



P H g C p\ T b,n~ T mi\ 



H 





where the various fluid properties are those of air evaluated at the mean film temperature: 



(132) 



^ 




(133) 



The internal convective heat transfer coefficient is a function of internal glazing layer surface temperature, 
T bn , for the case of natural convection so it is necessary to update the value of h cvmt as the solution of the 
glazed area heat transfer model proceeds. 



Each of the following Equations (134) to (139) pertains to one particular value, or range, of tilt angle, y. This 
categorization, as a function of / is based on the assumption that the internal environment is warmer than the 
internal glazing surface (i.e., ^jnt > ^b n)- ^ * ne reverse is true (Tjnt < 7" b n ) it is necessary to seek the 
appropriate correlation on the basis of the complement of the tilt angle, 180° - y, instead of y, and to then 
substitute 180° - y instead of /when the calculation is carried out. 

a) Windows inclined from 0° to 15° (0°< y< 15°) 




Nu- mi =0,-\3Ra H V3 



b) Windows inclined from 15° to 90° (15° < y°^ 90°) 
arrSinr ;RaTT^Ra r u 




1/4 



,56(_/?a cv sin/) ;Ra H >Ra 



c) Windows i 



o,72 7 y/ 5 

— ; /in degrees 
sin/ 

m 90° to 179° (90° < /< 179°) 



(134) 

(135) 
(136) 

(137) 



s inclined fror 
N M i nt =0,56(7?a i/ sin/)"' t ;10 o < Ra H sin/ <10 ' 
d) Windows inclined from 179° to 180° (179° < /< 180°) 



NOTE 



Nu- mt = 0,58Ra]i 5 ; Ra H ^ 10 11 



For more information on the convective heat transfer coefficient of still air see [27]. 



(138) 



(139) 
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8.3.2.3 Forced convection (any tilt) 

The following relation is to be used for the case of forced airflow on the internal side of a fenestration system. 
Vint=4 + 4F S (140) 

where V s is free stream velocity near the fenestration surfaces in metres per second. 
NOTE Equation (140) is taken from ISO 6946I 37 !. 

8.3.3 Convective heat transfer coefficient — external side 

8.3.3.1 General 

The convective heat transfer on the external side primarily occurs by forced convection. For those situations 
where natural convection occurs, see 8.3.3.5. The density of convective heat flow on the external boundary is 
defined as: 

<7cv,ex = "cv.ex (/s,ex ~* ex j (141) 

where T s ex is the temperature of any external fenestration surface (i.e., r f or the temperature of the external 
frame surface). 

8.3.3.2 Different applications 

There are two different applications that need to be considered for fenestration outside convective heat 
transfer coefficient correlations: 

a) fenestration product comparisons (rating), see 8.3.3.3; 

b) real building (field situation) fenestration component annual energy analysis, see 8.3.3.4. 

8.3.3.3 Convective heat transfer coefficient correlation for product comparison or rating 

The following relation is to be used for forced convection on the external side of a fenestration system: 

Vex=4 + 4F S (142) 

8.3.3.4 Real building fenestration component annual energy analysis 

Windows are major factors in overall energy performance of the real buildings. In order to be able to estimate 
the contribution to a building's thermal balance, it is necessary to define heat transfer characteristics of glazing 
as a function of different climate variables. Forced convective heat transfer on the external side of the building 
occurs between the surface of the building and the surrounding air at a rate that is determined by several 
factors. These include the temperature difference between the surface and the air, the speed and direction of 
any air movement (wind) over the building, and the shape and roughness of the building surface. Since those 
factors are highly variable, an exact mathematical analysis of the external surface convective heat transfer is 
not possible at this time due to the difficulty in defining the external surface and the highly variable wind 
conditions. 

NOTE For more information on this procedure, see Reference [28]. The correlations are based on the building 

experimental work given in Reference [29]. 
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Calculation procedure: 



^cv,ex=4,7 + 7,6F s 



(143) 



where 



if the surface is windward: 
F s = 0,25 V; V>2 m/s 



F s = 0,5; F^2m/s 



(144) 
(145) 



where V is the wind velocity measured at a height of 10 m above ground level and V s is free stream velocity 
near the fenestration surfaces; 



if the surface is leeward: 



V s = 0,3 + 0,05 V 



(146) 



To determine whether the surface is windward or leeward, calculate the wind direction, / w , relative to the wall 
surface (see Figure 26): 



rw = raz+ 180 °-rwN 

If \yj >180°,then rw = 360°- \yj 

If - 45° < | / w | < 45°, the surface is windward, otherwise the surface is leeward. 



(147) 




Key 

1 building 

n wall normal direction 

TwN w ' nc ' direction (angle measured clockwise from north) 

/ az wall azimuth (positive degrees westward from south and negative eastward) 

N north 
S south 

Figure 26 — Determination of wind direction and wall azimuth 
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8.3.3.5 Heat transfer by natural convection 



The natural convection heat transfer coefficient for the external side, h cvex , is determined in terms of the 
Nusselt number, Nu. 

A 



where X is the thermal conductivity of air. 

Nu is calculated as a function of the corresponding Rayleigh number based on the height, H, of the glazing, 
Ra H . 

P H gC p \T sex - T e J 
Ra H = 7p—, ~ ( 149 ) 

where the various fluid properties are those of air evaluated at the mean film temperature: 

Vf =7 ex+^s,ex-7ex) (150) 

Correlations to quantify the external side convective heat transfer coefficient are identical to the ones for 
internal side and are presented in 8.3.2.1. The tilt angle y needs to be replaced by its complement angle 
180° -/. 

NOTE For more information, see [27]. 

8.4 Longwave radiation heat transfer 

8.4.1 Mean radiant temperature 

External mean radiant temperature will depend on the application, whether it is for field conditions or for 
product rating and comparison (i.e., controlled laboratory conditions). For field conditions, external irradiance 
can be defined through the use of external mean radiant temperature, T m ex : 

-^ex = cr -'rm,ex 051) 

It is assumed that external fenestration surfaces are irradiated by the external surfaces and the sky vault 
consists of two areas, one cloudy and the other clear. The cloudy portion of the sky is treated as large 
enclosure surfaces existing at the external air temperature. The mean radiant external temperature can then 
be defined as: 



rm.out 



(^gd +( 1 ~/clr)' F sky) cr7, ex + fc\r F sy J sty 



(152) 



where F„ d and F sk are view factors from the external surfaces of the fenestration system to the ground 
(i.e., the area below the horizon) and sky, respectively. The factor/ dr is the fraction of the sky that is clear. 



^gd=1-^sky (153) 

77sk ^1 ± cos^ (154) 



If the radiosity of the clear sky (J sky ) is known, it can be used directly in Equation (152). 
NOTE Alternatively, if actual sky data are unavailable, the model from Reference [36] can be used. 
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-'sky — ^sky^ex U^v 

^sky=^ (156) 

R sky =5,31x10" 13 r 6 (157) 

Internal irradiance is defined as: 

^int = ^rm.int (158) 

where T m int is determined from temperatures and shape factors of surrounding internal surfaces. 

It is often assumed that internal fenestration surfaces are irradiated only by the internal room surfaces, which 
are treated as a large enclosure existing at the internal air temperature. Internal irradiance then becomes: 

£int=^int (159) 

The procedure outlined in this clause can be adapted to account for conditions that exist in a hot box test 
apparatus by determining the radiosities of the surfaces to which the window is exposed and the 
corresponding shape factors. 

8.4.2 Detailed radiation heat transfer calculation 

8.4.2.1 General 

Fenestration systems whose ratio of total to projected boundary surface area on the internal/external side is 
greater than 1,25 are called non-planar fenestration systems. For these systems, individual fenestration 
surfaces (i.e., frame and glazing surfaces) are self-radiating and the assumption of large black body enclosure 
radiating at each fenestration surface with the view factor equal to 1,0 is invalid. This analysis may also be 
used for certain fenestration components, such as frame cavities and ventilated cavities and grooves. 

The net radiation heat transfer on fenestration boundaries, q r , of non-planar products shall be calculated using 
procedure outlined in 8.4.2.2 or the alternative procedure given in 8.4.2.3. 

8.4.2.2 Two-dimensional element to element view factor based radiation heat transfer calculation 

The emissivity of both internal and external environments is set to unity. 

The net radiation heat transfer at any surface "/" is the difference between emitted radiation and absorbed 
portion of incident radiation. The temperatures of the surfaces do not appreciably differ so, using Kirchhoffs 
law: 

q ri =s i aT i -£iE t (160) 

where E { is irradiance at surface i from all other surfaces. 

N 

E i = H F i-j J j ( 16i ) 

j 
and F t ■ is the view factor from surface i to surfacey. The radiosity of surface j, Jj, is given by: 

Jj=e j oT j A +Pj Ej (162) 
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Assuming all surfaces are grey: p- = 1 - y a .. Substituting p. and E- and using subscript / for convenience, 
Equation (162) becomes: 

N 

J t = Si aT 4 + (1 - e^Fi.jJj (163) 

7=1 

Equation (163) represents a system of TV linear algebraic equations for the TV unknown radiosities, J, which are 
determined from the solution of this system of equations. The system of Equation (163) when expressed in 
matrix form becomes: 

[C]{J}={D} (164) 

where 

= _jl_A LT_1^. (165) 

s t 

D { =gT 4 (166) 

Tj in Equation (166) is a known temperature from previous iteration k, (i.e., T t \ k ). For the first iteration, the 
values for T { are initial guesses. 



Temperatures are calculated from the solution to the conduction problem given by Equation (72), while net 
radiation heat flow rate [see Equation (160)] is calculated using J t values from Equation (163) and linearized 
term T t 4 , by using the first two terms of its Taylor series expansion about T i 



\k 



T 4 ^( Ti \ k f Ti \ M ~^i\ k ) 4 (167) 

This procedure is repeated until the following condition is satisfied: 

\k 



rp\k+'\ 



r |*+1 



<to/ (168) 



where tol is solution tolerance, whose value is typically less than 10~ 3 . | | denotes the norm or root mean 
square value of the temperature vector. 

View factors F u - can be calculated using Hottel's cross-string rule. If the view between two radiating surfaces 
is obstructed by a third surface, the effect of this obstruction shall be included. 

NOTE For more information on the norm or root mean square values of the temperature vector, see [30]. Also, for 

more information on the cross-string rule, see [26]. 

8.4.2.3 Simplified three-dimensional radiation method 

The thermal transmittance of non-planer products must be reduced because of the self-viewing nature of the 
window. The alternative method presented here can be used instead of the multi-element method prescribed 
in 8.4.2.2. 

The emissivity of the internal surface is reduced by the factor given in Equation (169). The ratio of internal- 
side surface area to internal-side-opening area is denoted as A s /A p . The opening area, A p , will be similar to, 
but marginally less than, the projected area of the product, A t . 

NOTE More detail about this method is provided in Reference [31]. 
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F r = 7 1 (169) 



1 + ^g 



A, 



In the analysis of individual window glazing segments, the radiant exchange from the internal glazing surface 
can most readily be reduced by substituting a reduced emissivity for the internal surface, £- red , in place of the 
true surface emissivity, s b n . 

^red=^r^b,« ( 170 ) 

Similarly, in the analysis of frame and sash sections, the emissivity of the internal surfaces shall also be 
reduced using the factory^. 

8.4.3 Simplified radiation heat transfer calculation 

8.4.3.1 Internal surfaces 

All internal surfaces are denoted by subscript s int , including frame surfaces. Although it is customary to 
consider radiation heat transfer on glazing boundary surfaces in terms of radiosities, as shown in 5.3.1, the 
following equation can be used for simplified radiation heat transfer calculations on both glazing and frame 
surfaces: 

9r,int = Vint(,7s,int -7 Vm,int) (171) 

where 

£ s, int " 's, int _ -'rm,int 

Vint= =} " L (172) 

J s,int -'rmjnt 

8.4.3.2 External surfaces 

All external surfaces are denoted here by subscript sex , including frame surfaces. As with internal surfaces, 
external glazing boundary surfaces are typically dealt with in terms of radiosities, but the following equation 
can be used for simplified radiation heat transfer calculations on both glazing and frame surfaces. 

<7r,ex = "r,ex (/ s,ex - ^rm,ex j V'^) 

where 

£ s, ex " 's, ex _ -'rm,ex 

Vex= ^ " L (174) 

1 s,ex 1 rm,ex 

8.5 Combined convective and radiative heat transfer 

q = h(T s -T n ) (175) 

where 

h = h r +h c (176) 

T s is the surface temperature; 

T n is the environmental temperature. 
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For a boundary surface whose geometry is approximated using rules described in 6.3.1, 
adjustment to the combined surface heat transfer coefficient shall be applied: 



"adjusted 



'real 



^approximated 



-x h 




the following 
(177) 



8.6 Prescribed density of heat flow rate 

The frame/wall interface shall be treated as adiabatic. Also see ISO 10077-2. 



To calculate total solar energy transmittance, r S |, of frames, g f in analogy to Equation (14) 
boundary condition shall be used: 



9int =aI s 



where a is the absorption of frame surfaces 
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Annex A 

(informative) 

Solution technique for the multi-layer solar optical model 



A window with n glazing together with the external (i = 0) and internal (i = n + 1 ) spaces form an n + 2 element 
array. The glazing system optical analysis can be carried out by considering the spectral fluxes of radiant 
energy flowing between the / - 1th and z'th glazings, //"(A) and If (A). The + and - superscripts denote 
radiation flowing toward the external and internal side, respectively, as shown in Figure A.1. 






r 

i 


>7 + i 












'"ft,/ 

— ^ — 


r b./ 

— (A) 

7; 



n + 1 



/ + 

n + 1 



/-1 



/+ 1 



Figure A.1 — Analysis of solar flow rate in multi-layer glazing system 

Equations A.1 and A.2 are applied while setting the reflectance and transmittance of the conditioned space to zero, 
r fti1 (A) = and r„ +1 (A) = 0. 



/ r + (A) = r i (A)/ + 1 (A) + r fu (A)/r(A) 



1 to n + 1 



(A.1) 



/"(A) = r,_ 1 (A)/- 1 (A)+r b|i _ 1 (A)/, + (A) 



/ = 2 to n + 1 



(A.2) 



It can be shown that the ratio of 7,_-|(/l)to /,- (A) is given by: 






T?(A)r M (A) 






(A3) 



and the ratio of 7 r _ 1 (/l)to I, (A) is given by: 






■tM): 



tM) 



1-r b (A)r /+1 (/l) 



(A.4) 



NOTE Reference [32] gives the development of Equations (A.3) and (A.4) by the net radiation method. Reference [5] 

uses ray tracing to reproduce this derivation. 

Now all values of I[(A) and 7 r + (/l) can be found using the following steps. First, use Equation (A.3) as a 
recursion relationship to calculate all values of r t {X) by working from ^ + i(A) = r ft „ +1 (A) = to r n _ x {X). 
Second, use Equation (A.4) to obtain r,-(A) from i=\ to i = n. The rest of the solution follows by first 
calculating the density of heat flow rate reflected from the glazing system to the external side [i^(A) = r^{A)i^(A)] 
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and then marching toward the internal side from / = 2 to / = 1 calculating the remaining density of heat flow 
rates [/f(A) = r / _ 1 (A)/ < i 1 (A) and i+(X) = n (X)if(X)]. 

Finally, the desired portions of incident radiation at each of the glazing layers are given by: 

A (,l)- A r W-^ + W + ^+iW-*r + iW (A.5) 

and the portion transmitted to the conditioned space is: 

w Jn^W (A.6) 

7f(A) 

7f (A) can be set equal to unity for the purpose of solving these equations. 

The approach used here can also be applied to the analysis of shading devices (i.e., Clause 7) although the 
resulting equations will be more complex due to the conversion of beam solar radiation to diffuse solar 
radiation. 



©ESA 61 



ESISO15099:2012(E) 



Annex B 

(normative) 

Thermophysical fill gas property values 



The following tables present linear equation coefficients with which the thermal conductivity, dynamic viscosity 
and specific heat capacity at constant pressure can be calculated for four (air, argon, krypton, xenon) glazing 
cavity fill gases. The least square linear curve fit equations were derived from thermophysical property data 
given in Reference [33]. 

The heat transfer calculations are based on the assumption that the fill gas is not an emitting/absorbing gas. 
Since SF 6 violates this condition it has not been presented in these tables. 





Table B.1 - 


— Thermal conductivity 




Gas 


Coefficient a 

W/(mK) 


Coefficient /> 

W/(m-K 2 ) 


latO°C 


/lat10°C 


W/(m-K) 


W/(m-K) 


Air 


2,873x10- 3 


7,760 x10- 5 


0,024 1 


0,024 8 


Argon 


2,285x10- 3 


5,149 x10- 5 


0,016 3 


0,016 9 


Krypton 


9,443 x10- 4 


2,826 x10- 5 


0,008 7 


0,008 9 


Xenon 


4,538 x10- 4 


1,723 x10- 5 


0,005 2 


0,005 3 


where X= a + b-T(K), in 


W/(mK). 









Table B.2 


— Dynamic viscosity 




Gas 


Coefficient a 

N-s/m 2 


Coefficient b 

Ns/(m 2 K) 


juat0°C 


juat10°C 


Pas 


Pas 


Air 


3,723 x10- 6 


4,94 x 10- 8 


1,722x10- 5 


1,771 x10-5 


Argon 


3,379 x10- 6 


6,451 x10" 8 


2,100 x10- 5 


2,165x10-5 


Krypton 


2,213 x10- 6 


7,777x10- 8 


2,346x10- 5 


2,423x10-5 


Xenon 


1,069 x10- 6 


7,414x10- 8 


2,132x10-5 


2,206x10-5 


where fj = a + b-T{K), 


n Pas. 








NOTE 1 Pas = 1 


kg/(ms) = 1 N-s/m 2 . 
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Table B.3 — Specific heat capacity at constant pressure 



Gas 


Coefficient a 

J/(kgK) 


Coefficient b 

J/(kgK 2 ) 


c p at °C 
J/(kgK) 


c p at10°C 
J/(kgK) 


Air 


1 002,737 


1,232 4 x 10- 2 


1 006,103 4 


1 006,226 5 


Argon 


521,928 5 





521,928 5 


521,928 5 


Krypton 


248,090 7 





248,090 7 


248,090 7 


Xenon 


158,339 7 





158,339 7 


158,339 7 


where c = a + b-T(K) in J/(kgK). 



Table B.4 — Molecular masses 



Gas 


Mass 

kg/kmol 


Air 


28,97 


Argon 


39,948 


Krypton 


83,80 


Xenon 


131,30 
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Annex C 

(informative) 

Examples of calculated values for optical properties of slat type 

of shading devices 



C.1 General 

This annex provides selected examples of calculated optical properties of a slat type of shading device as 
function of the slat properties and geometry, for different angles of incident solar radiation, applying the 
procedure given in 7.3. 



C.2 Specifications 

See Figures C.1 and C.2. 
Geometry: 

slat distance: 12 mm 
slat length: 16 mm 
slat angle: 45° and/or 80° 

Solar incidence angle: 

0° and 60° from horizontal 
Slat types: 

opaque or translucent 

white, pastel and/or dark coloured (outside and/or inside surfaces) 
The optical properties of the materials of the different types of slats are given in Table C.1 . 

Table C.1 — Optical properties of the slat materials 



Slat 


Long wave (IR) 
transmittance 

T 


Long wave (IR) 
emissivity 


Solar transmittance 


Solar reflectance 

r 


Opaque, white 


0,0 


0,90 


0,0 


0,70 


Opaque, pastel 


0,0 


0,90 


0,0 


0,55 


Opaque, dark 


0,0 


0,90 


0,0 


0,40 


Translucent, white 


0,40 


0,55 


0,40 


0,50 
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Dimensions in millimetres 






Key 

1 slat 

2 solar angle 




Figure C.1 — Slat geometry and sun position 




vV 




Figure C.2 — Slats angle 45° and 80° 
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C.3 Results 

The equations in 7.3 are used to calculate the optical properties of the shading device. 

The optical properties are illustrated in Figures C.3 and C.4 and the results of the calculations are presented 
in Table C.2. 

Examples of direct-direct transmittance r dirdir 
Radiation from outdoor side Radiation from indoor side 

Opaque slats Translucent slats Opaque slats Translucent slats 




















<? 























T dir, dir 












































1 







2 




Key 

1 outdoor 

2 indoor 

Figure C.3 — Illustration of the direct-direct transmittance of the shading device 

Examples of direct-diffuse transmittance r dir dif and Examples of diffuse-diffuse transmittance r dif dif and 
reflectance r dir dif reflectance r dif dif 

Radiation from outside Radiation from inside Radiation from outside Radiation from inside 




Key 

1 outdoor 

2 indoor 



Figure C.4 — Illustration of the direct-diffuse and diffuse-diffuse transmittance 
and reflectance of the shading device 
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Note that the optical properties for radiation from the inside are in general not the same as the optical 
properties for radiation from the outside. The optical properties for radiation from the inside are needed 
because of (multiple) reflections at the successive layers in the transparent system. 



Table C.2 — Selected examples of calculated optical properties of 


a slat type of shading 


device 


Product 


Description 


A45 


Venetian blinds with opaque white slats, slat angle 45° 


B45 


Venetian blinds with opaque pastel colour slats, slats angle 45° 


C45 


Venetian blinds with opaque light/dark slats, slats angle 45° 


C80 


Venetian blinds with opaque light/dark slats, slats angle 80° 


D45 


Venetian blinds with translucent white slats, slats angle 45° 


Property 


A45 


A45 


B45 


B45 


C45 


C45 


C80 


C80 


D45 


D45 


Geometry: 




Distance between slats 
(mm) 


12 




12 




12 




12 




12 




Slat length (mm) 


16 




16 




16 




16 




16 




Slat angle, from 
horizontal (°) 


45 




45 




45 




80 




45 




Optical material 
properties: 




IR transmittance slat 


0,0 




0,0 




0,0 




0,0 




0,40 




IR emissivity slat 
external side 


0,90 




0,90 




0,90 




0,90 




0,55 




IR emissivity slat 
internal side 


0,90 




0,90 




0,90 




0,90 




0,55 




Solar transmittance slat 


0,0 




0,0 




0,0 




0,0 




0,40 




Solar reflectance slat 
external side 


0,70 




0,55 




0,70 




0,70 




0,50 




Solar reflectance slat 
internal side 


0,70 




0,55 




0,40 




0,40 




0,50 




Solar incidence angle 
from horizontal (°) 





60 





60 





60 





60 





60 


Results: 




Solar transmittance 
external side dir, dir 

» dir.dir'ex 


0,057 


0,0 


0,057 


0,000 


0,057 


0,000 


0,000 


0,000 


0,057 


0,000 


Solar transmittance 
internal side dir, dir 

( r dir,dir)int 


0,057 


0,310 


0,057 


0,310 


0,057 


0,310 


0,000 


0,088 


0,057 


0,310 


Solar transmittance 
external side dir, dif 

( r dir,dif)ex 


0,141 


0,073 


0,090 


0,047 


0,096 


0,051 


0,012 


0,005 


0,373 


0,277 


Solar transmittance 
internal side dir, dif 

( r dir,dif)int 


0,141 


0,288 


0,090 


0,216 


0,076 


0,271 


0,011 


0,027 


0,373 


0,306 
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Table C.2 (continued) 



Property 


A45 


A45 


B45 


B45 


C45 


C45 


C80 


C80 


D45 


D45 


Solar reflectance 
external side dir, dif 

( r dir,dif)ex 


0,394 


0,558 


0,295 


0,430 


0,371 


0,544 


0,622 


0,678 


0,418 


0,567 


Solar reflectance 
internal side dir, dif 

( r dir,dif)int 


0,394 


0,103 


0,295 


0,066 


0,216 


0,070 


0,356 


0,273 


0,418 


0,273 


Solar transmittance 
external side dif, dif 

( r dif,dif)ex 


0,332 


0,294 


0,291 


0,038 


0,495 


Solar transmittance 
internal side dif, dif 

( r dif,dif)int 


0,332 


0,294 


0,291 


0,038 


0,495 


Solar reflectance 
external side dif, dif 

(Aiif,dif)ex 


0,345 


0,260 


0,323 


0,604 


0,380 


Solar reflectance 
internal side dif, dif 

( r dif,dif)int 


0,345 


0,260 


0,193 


0,345 


0,380 


IR transmittance 
external side (zj R ) ex 


0,227 


0,227 


0,227 


0,0245 


0,385 


IR transmittance 
internal side (zj R ) int 


0,227 


0,227 


0,227 


0,0245 


0,385 


IR emissivity 
external side (s) e>dernal 


0,729 


0,729 


0,729 


0,89 


0,536 


IR emissivity 
internal side (s) mt 


0,729 


0,729 


0,729 


0,89 


0,536 
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